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High-Temperature Thermodynamic Functions for 
Zirconium and Unsaturated Zirconium Hydrides 

Thomas B. Douglas 
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Giving greatest weight to the experimentally measured highest decomposition pressures 
and the enthalpies in one-phase fields, thermodynamically interconsistent integral and 
differential enthalpies (heat contents), heat capacities, entropies, and Gibbs free energies are 
derived for the crystalline one- and two-phase fields of the zirconium-hydrogen system for all 
stoichiometric compositions from Zr to ZrHi.25 and over the temperature range 298.15 to 
1,200 °K. These properties are derived in analytical form, and in most cases are represented 
by numerical equations, with tabulation for zirconium and H/Zr atom ratios of 0.25, 0.50, 
0.57, 0.75, 1.00, and 1.25. Most of the unique phase-field boundaries which are consistent 
with the derived properties are located and are compared with those previously reported. 
In the Zr-H system the enthalpies are shown to relate certain properties at different com- 
positions as well as at different temperatures. Some of the various data show good inter- 
consistency, while others reveal discrepancies which are discussed critically. 
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1. Introduction 

Among transition-metal hydrides, the zirconium- 
hydrogen system has been studied very extensively 
in recent years because of the practical importance 
of and theoretical interest in the system. Zirco- 
nium hydrides are potentially valuable thermal- 
neutron moderators which in some cases can reach 
elevated temperatures without attaining excessive 
pressures. The wide ranges of composition, in- 
terrupted by lattice changes and marked miscibility 
gaps, offer a fruitful system for studying the structure 
of interstitial solid hydrides. 

Experimental data on the thermodynamic and 
other physical properties of the system have been 



published by a number of investigators. Some of 
these have measured equilibrium pressures of hydro- 
gen in certain ranges of temperature and composition 
II, 2, 3, 4, 5, 6, 7, 8]. 1 The rates of diffusion of 
hydrogen in the metal have been studied [9, 7]. 
X-ray |10| and neutron-diffraction [11] studies have 
given information on the zirconium- and hydrogen- 
atom spacings, respectively. Measurements have 
been made of thermal conductivity [12] and thermal 
expansion [12, 13, 14]. The enthalpies or heat 
en pa cities have been measured for the "saturated" 
composition ZrH 2 at low temperatures [15], and for 
several "unsaturated" compositions at high tem- 
peratures [16, 12]. Many of these data have served 
to outline the phase diagram of the system by a 
variety of methods, as the author showed in a 
previous paper [17]. Attention may also be called 
to some of the publications which have been par- 
ticularly directed toward reviewing these properties, 
correlating them thermodynamically, and discussing 
their bearing on an interpretation of the structures of 
these phases [4, 15, 17, 18, 19]. 

In this paper is presented a consistent set of the 
common integral and partial thermodynamic func- 
tions of the zirconium-hydrogen system based on 
some of the aforementioned data. These functions 
were derived from the data without recourse to 
extra-thermodynamic considerations except, inci- 
dentally, such general and well-established ones as the 
laws of dilute solutions. The composition range 
covered is from Zr to ZrH 12 5, and the temperature 
range from 298.15 to 1,200 °K. (Though thermody- 
namic data are available for higher temperatures [1] 
and higher hydrogen contents [2,3,5,8,10], these were 
omitted from consideration because of their generally 
inferior precision or lack of agreement.) The co- 
existence of phases in equilibrium is in a sense acci- 
dental, and the theoretical interest in the system 



1 Figures in brackets indicate the literature references at the end of this paper. 
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therefore centers mainly in the properties of indivi- 
dual phases. But from a practical standpoint the 
equilibrium reproportionation of phases with temper- 
ature change of a given overall composition occurs 
fairly rapidly in some cases and is of interest; for 
this reason properties have been evaluated for two- 
phase as well as for one-phase regions. As will be 
pointed out, the available data show excellent inter- 
consistency in some ranges of composition and tem- 
perature. In other regions considerable discrepancies 
were found which introduce considerable uncertainty 
into the adopted thermodynamic functions, partic- 
ularly in two-phase regions. It will be evident that 
the requirement of thermodynamic consistency for a 
multicomponent system imposes an often unsus- 
pected multitude of restrictions on the relations of 
the properties for neighboring compositions and tem- 
peratures, particularly near the center of the region 
under consideration in the phase diagram. This is 
true even for a two-component system which, like 
the present one, has a phase diagram which is only 
moderate!}^ complex. 

The opinion has been expressed [6] that the purity 
of the zirconium used in the various experimental 
studies on the zirconium-hydrogen system was not 
sufficiently high for the results of a full thermody- 
namic analysis to be of any great significance. The 
present paper will illustrate that such an analysis 
can be rather laborious, and that, in the present case, 
it reveals numerous inconsistencies in the data which 



lend some weight to that opinion. Nevertheless, the 
author feels that only by a thermodynamic correla- 
tion of this type can the interconsistencies of the 
various data be tested, a better idea of their reli- 
ability be gained, additional properties of the system 
be derived, and the need for additional experimental 
work in certain areas be suggested. However, to the 
extent that the data are inconsistent, an adopted 
correlation will necessarily be somewhat arbitrary, 
particularly since the theoretical number of different 
ways to weight the numerous interrelated data is 
found to be quite large. In another sense, one can 
be arbitrary in deciding which borderline theoretical 
relationships are sufficiently valid to impose on the 
correlation. 

The phase diagram of the zirconium-hydrogen 
system in the region treated in this paper is shown in 
figure 1. The symbols a, /3, and d 2 refer to single- 
phase regions where the zirconium lattice is close- 
packed hexagonal, body-centered cubic, and face- 
centered tetragonal/ respectively. The points 
shown on the diagram are those on phase-field 
boundaries as found by various investigators, and 
the solid boundaries are those to which the tliermo- 



2 Some authors have used "y" in place of "<5" in this system, "y," however, 
has come to be used [3, 12] to designate a type of metastable phase in the system 
not discussed in the present paper. 

3 Face-centered cubic may here be considered a special case of face-centered 
tetragonal. By a simple rotation of crystal axes, the 6 face-centered tetragonal 
lattice may be transformed into body-centered tetragonal (with a different unit 
cell) , for comparison with the /3 body-centered cubic lattice. 
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Figure 1. Phase diagram (partial) of the zirconium-hydrogen system. 



(Observed points on phase boundaries are shown as follows. From equilibrium hydrogen pressures: Q Hall, Martin, and 
Rees [2]; □ Gulbransen and Andrew [41; G2 Edwards, Levesque, and Cubicciotti [5]; O Ells and McQuillan [6]; X Mallett 
and Albrecht [7]; <g> LaGrange, Dykstra, Dixon, and Merten [8] From hydrogen diffusion: V Schwartz and Mallett [9]. 
From X-ray diffraction: AVaughan and Bridge [10]. From thermal expansion: I] Beck [12]; H Espagno, Azou, and Bastien 
[14]. From enthalpy: O Douglas [17]. The interconsistent thermodynamic functions formulated in this paper correspond 
to the solid-curve boundaries, and are tabulated along the dotted lines.) 
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dynamic properties adopted in t his paper correspond. 
The majority of these properties are presented in 
analytical form for the continuous ranges of tem- 
perature and composition, and are tabulated along 
the doited vertical lines shown in the figure. 

2. Thermodynamic Preliminaries 

2.1. Conventions 

The unit of energy used throughout is the defined 
thermochemical calorie (=4.1840 j), and all tem- 
peratures are in deg K. All extensive properties of 
a zirconium hydride are evaluated for "one mole/' 
defined as containing a total of 1 g-atom of Zr and x 
g-atoms of H, 4 even when one mole so defined is com- 
posed of two or more phases in equilibrium. The 
components are taken to be alpha-zirconium (at all 
temperatures, to preserve continuity into the small 
64° temperature interval covered where beta- 
zirconium is more stable) and "normal" molecular 
hydrogen in the usual standard state of the ideal 
gas at 1 atm. 5 

The symbols in common usage have been adopted 
so far as possible. Those not specifically defined 
are : 

di\ Thermodynamic activity of com- 
ponent i. (a 2 is taken equal to the 
equilibrium pressure of H 2 (g) in 
atm.) 
h, c, b+c: Generalized notation for particu- 
lar one- or two-phase fields. (In 
the absence of specific designation, 
the phase field is implied as any 
one, or the particular one under 
discussion.) 
C v \ Heat capacity at constant pres- 
sure and constant overall compo- 
sition (cal mole -1 deg K" 1 ) . 
/: "Of formation" (usually from the 4 
components in their standard 
states) . 
G: Absolute Gibbs free energy (cal 

mole -1 ). 
H: Absolute enthalpy (heat content) 

(cal mole -1 ). 
H'\ "Intrinsic enthalpy" (AH of the 
reaction of eq (30)). 
H 2 (g) : Hydrogen gas. 
i (subscript): A particular component (either). 
In: Natural logarithm. 
log: Logarithm to the base 10. 
R : The molar gas constant ( = 1 .98725 

cal mole -1 deg K _1 ). 
S: Absolute entropv (cal mole" 1 deg 

K- 1 ). 
T: Absolute temperature (deg K, Int. 
Temp. Scale of 1948). 



' This convention si 
it is probable that liy< 

(willi vacancies) Insto 

'• By this con \ mi io 

moles of the compom 

and atomic IT could o 



more logical t han one based on mole fractions, since 
(as atoms) enters the zirconium lattices interstitially 
substitutional^ 7 . 

mole of ZrHj thus contains one mole of Zr and \ x 
The conversion to standard states involving /S-Zr 
e readily be made if desired. 



x: Atoms of H per atom of Zr (in the 
overall composition when two 
phases are present). 
x b and x c : Define the compositions of the two 
coexistent phases ZrH^ and ZrH^ 
which lie on the phase-field bound- 
aries of the (b + c) field. 
Y: Generalized notation for one of H, 
_ C PJ S, or G. 

Y t : Partial niolal Y with respect to 
component i. 
a, P, 8: Types of single phases of ZrH x . 
(See fig. 1.) 
(a-f-/?), etc.: Two-phase field, or an overall com- 
position lying therein. 
a/(a-\-ff) y etc.: Boundary between the a and 
(a + ff) fields, etc, 

1 (subscript): Tlie component zirconium, Zr. 

2 (subscript): The component molecular hydro- 

gen, Ho. 
A numerical subscript or superscript other than 1 or 
2: When ^2, a particular value of x; when >2, 
a particular value of T. (The absence of such 
subscripts implies general values of x and T, or a 
particular set under discussion.) 
or ° (superscript): In the standard state (see 
text); deg temperature. 

2.2. General Thermodynamic Relations 

Several well-known thermodynamic relations fre- 
quently used in subsequent derivations or calcula- 
tions are assembled here. These have been adapted 
to the specialized definition of one mole of Zr, H 2 , 
and ZrHz stated above. 

Total properties are often conveniently derived 
using 

C p =[d(H -H° 298 , 15 )/dT) x] (1) 

(dS°/dT) x =C;/T; (2) 

-(G o -H% 8 . 15 )/T=S -(H°-ir 29S . 15 )/T. (3) 

(Equation (3) has been used throughout this paper 
to obtain the free-energy function.) A relation 
similar to eq (3) applies to the relative partial molal 
properties: 

(6 ( - G°)/T= (Ht-HV/T- (S t -S°) . (4) 

The total and partial properties are related by 



and 



Y°=Y l +izY t 



Y 2 =2{bY°lbx) T , 



(5) 



(«) 



integration of the latter between two compositions 
x f and x" giving 



fe A T 111 & 



TO TrO 



Y° „— I 



1 C £ " — 



(7) 
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From eqs (5) and (6) 

F 1 =r°-x(dF°/^) T . (8) 

The Gibbs-Duhem equation takes the form 

[d(F 1 -FJ)/Sa;] r +Ma(F 2 -F°)/^] r =0. (9) 
By the definition of activity 

(G-G°^IT=Rh\a i . (10) 

The variation of the partial molal free energy with 
temperature is given by 

[b(^=^)/dr],= -(F 4 -H?)/^. (id 

From eq (5) follows the standard heat, entropy, 
and free energy of formation at any temperature. 
Since 

AYf=Y°-(Y° 1 +ixY° 2 ), 



we have 



AHf= fa-Hi) + hx(H z -Hl) ; 

AS f=(S 1 -S° 1 ) + kx(S 2 -S? 2 ); 
AGf=(G-G° i ) + MG 2 -Gl). 



(12) 

(13) 
(14) 
(15) 



The molal enthalpy, heat capacity, entropy, or 
free energy in a two-phase field is the sum of those 
of the component (or "terminal") phases: 

Yl+°=[(x-x)l{x-x b )]Y»+[{x-x b )l(x-x b )]Y°. 

(16) 

In the present case it is particularly simple to relate 
two common alternative conditions determining 
this separation into the two terminal phases. One 
commonly stated obvious condition is that a straight 
line be tangent to the continuous free-energy iso- 
therm at x b and x c and lie lower than the isotherm 
between x b and x c . Hence 



(dG!dx) x =(dG/dx) Xc , 
and hence from eq (6) 

(G 2 )x b : =(G2)x c =G2. 

From eqs (5) and (18) 



G x =(G 1 ) Xh +ix b G 2 



and 



tf* c =(#i)* c +W?2, 
whose difference gives 
[Gr-G Xh ]l[x-x,\=[m 



(17) 

(18) 

(19) 
(20) 



(6 , i)xJ/[* e -z»] + lG , 2. 



(2i; 



From the condition of tangency and eq (6) the first 
and third terms of eq (21) are equal, and hence 



(Gi)x b = (Gi) z 



(22) 



Equations (18) and (22) are the more commonly 
stated conditions for phase coexistence. The con- 
dition that the tangent line lie lower than the 
metastable free-energy curve between requires that 
at each point of tangency d 2 G/dx 2 be positive, and from 
equation (6) this is true if at each of these points 
dG 2 /dx is positive — i.e., if the activity of hydrogen 
in each phase increases isothermally with increasing 
hydrogen content. 

3. Basic Data Used 

The direct measurements of the equilibrium pres- 
sures of hydrogen for the system extend from 650 to 
1,373 °K and over the composition range from near 
x=0 to af>1.7. These data are sufficient to deter- 
mine a set of the common energetic thermodynamic 
properties of the system in this range. However, 
these sets of data overlap considerably and with 
varying degrees of disagreement. In addition, those 
of the precise enthalpy data which may be considered 
valid, while they alone neither establish activities 
nor their temperature derivatives, do impose rela- 
tions among the activity values at a given composi- 
tion and different temperatures and also, as will be 
presently shown, among the temperature derivatives 
of activity for different compositions. 

Enthalpy measurements relative to T=273 were 
made by Douglas and Victor [16] (usually at tem- 
peratures spaced at 50- or 100-deg intervals) up to 
T=l,173 for the compositions x=0, 0.324, 0.556, 
0.701, 0.999, and 1.071. These measurements were 
all made by the drop method, and in all cases in- 
volved fairly rapid cooling in the low-temperature 
two-phase (a +5) field, where composition equilib- 
rium corresponds to phase reproportionation (see 
fig. 1). However, evidence will be shown in section 
5 that phase equilibrium was approximately main- 
tained in this field in these measurements. It would 
also be expected that composition equilibrium was 
attained for temperature intervals beginning in the 
one-phase field involved (the field), and this was 
indicated by the high precision, regularity, and free- 
dom from hysteresis found for most of these 
particular measurements (see sec. 6). 

In view of the above facts, it was decided to 
construct the adopted consistent set of thermo- 
dynamic functions for the system largely by giving 
greatest weight to the hydrogen-activity data at the 
higher temperatures and to use the enthalpy data 
considered valid in order to extend the properties so 
derived to other temperatures and other composi- 
tions. The measurement of the highest equilibrium 
pressures is ordinarily subject to the smallest per- 
centage errors, and this is often true with respect to 
the effects of small amounts of interfering impurities 
in the samples. 

The hydrogen-activity data so chosen as a basis 
were those reported graphically by LaGrange, Dyk- 
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stra, Dixon, and Merten [8] at T= 1,073.2 and 13 
compositions (a=0.05 to 1), and (hose reported in 
tabular form by Edwards, Levesque, and Cubicciotti 
[5] for values of x of 0.650, 0.733, 0.815, and 0.894 
and at 8 temperatures over the range 2"= 1,023 to 
1,173. (These points all lie in the field.) The 
values of the latter authors were first all reduced to 
7=1,073.2 by using their smoothed temperature co- 
efficients in the usual way. A reasonably good fit 
to these data was found to be given by 



J 2 (1073. 2)- 



=0.06084x 136 +0.11251x 8 



(23) 



or 



[(«- 



■ A8)/rW 2=6.22313 log x 

+4.5758 log (l + 1.8493x 6 - 97 )- 



-5.56334. (24) 



The agreement with the values fit (which vary from 
0.002 to 0.16 atm) is shown in table 1; over the 

Table 1. Deviation from eg {23) {or (IB 4)) of observed hydro- 
gen activities of /3-ZrII, at 1,073.2 °K 



Reference 


X 


a 2 (obs)-a2(calc(l) x/ ](R) 
ct2(obs) 


[8 




0.05 


+43 


8 




. 11 


-12 


8 




. L9 





8 




. 27 


u 


8 




.37 


-1 


8 




.40 





8 




. 55 


+ 1 


[8 




.64 


() 


5 




. 650 


-4 


l- r ) 




.733 





8 




.76 





5 




.si 5 


+4 


[8 




.87 





[5 




.894 


+5 


[8 




. 92 


-4 


[8 




.96 





is 




1.00 


-8 



range #=0.19 to 1 it averages 2 percent, which is 
within the precision of the measurements. Calcu- 
lations of cycles of entropy and zirconium-activity 
which are presented in section 7 depend in part on 
the accuracy of eq (23) (or (24j), and show excellent 
consistency. 

It is appropriate to examine eq (23) in the light of 
the fact that several observers have found a pro- 
portionality, at constant temperature and only 
moderately high hydrogen concentrations, between 
the activity a 2 and the square of the hydrogen con- 
centration ("Sieverts' law"), which corresponds to 
an activity coefficient of atomic hydrogen independent 
of composition. Up to x=0.5 the last term in eq 
(23) contributes less than 2 percent of the total, 
so that t he power of xis not 2 but closer to 1. Sieverts 
and Roell [1] measured hydrogen-pressure isotherms 
for beta-zirconium at 1,073° and 1,373 °K, and 
though their data show inferior precision compared 
with that of the more recent measurements fitted 
above, their results show an approximate proportion- 
it lity of a 2 to x 2 at the higher temperature but not at 
the" lower. It is easy to show, by substituting from 
eq (10) into eq (11), replacing a 2 by a single term 
proportional to x n , and then differentiating partially 



with respect to x } that the power n is independent of 
T if and only if (H 2 —H° 2 ) does not vary with x. It 
will in fact be found that although (H 2 —H° 2 ) near 
1,073 °K is approximately independent of x at the 
higher values of ./\ it is far from being so for x<0.5. 
The sign and magnitudes of b(H 2 — H° 2 )/i>x in the latter 
region correspond to an approach to the a? 2 law as the 
temperature rises, but for lack of enough precise data, 
no serious effort was made to interpolate the pressure 
data over the range 1,073 to 1,373 °K. 

The manner of using the enthalpy data will next 
be considered. In the first place, it is apparent from 
figure 1 that in all these enthalpy measurements 
for the five hydride compositions the sample was 
cooled into the (a+d) field at 298 °K. In section 5 
the following value will be derived (from eq (90)): 



(Hr 8 -H°) 2QS =4. 



(25) 



This value is very small because at this temperature 
the composition range of_the a field is very small. 
We may now consider (H 2 +b — H° 2 ) 298 . Applying eq 
(6) to the absolute enthalpy at any temperature and 

also at 298 °K, and subtracting the two, one may 

write 



(ii 2 -ui)=2\b{ir-m« 
+(ii2 +8 



8) )/dx\ 



(26) 6 



The application of eq (6) to eq (16) shows that 
relative partial molal functions in a two-phase field 
at a given temperature (including (// 2 N — H° 2 ) 29H ) 
are independent of x. Consequently, if data at one 
temperature and composition are available for evalu- 
ating (H 2 +b — H 2 ) 298 from eq (26), then the same 
equation, with this value substituted, maybe used to 
evaluate (H 2 —H° 2 ) at other temperatures and com- 
positions where \d(II —H 29 g)/dx] T can be evalu- 
ated by suitable interpolation. 

Plots of the hydrogen-pressure data of Edwards et 
al. [5] for the four compositions from £=0.650 to 
^=0.894, used above, gave without trend almost 
identical values which averaged (H £ — H° 2 ) 1123 = 
—40,500. Ells and McQuillan [6] measured a? at one 
high concentration of hydrogen (z=0.66) from 
approximately 875 to 1,035 °K, leading to a value 
which when corrected to 1,123 °K by the enthalpy 
data gave (H§-H 2 ) 1123 = -37,500. A weighted 
mean was adopted: 



(Hi- 



^2) 1123,0.' 



-39,500. 



(27) 



For many calculations in this paper the relative 
enthalpy of molecular hydrogen gas has been com- 
puted from the empirical equation 



H o 2 -#°2(298)=6.660r 



+2.678(10- 4 )r 2 -12000r- 1 -1969, (28) 
which between 298 and 1,200 °K agrees within an 



i It as assumed in eq (26) and related eqs (30-36, 39, 40) that the values of x lie 
in the («+5) field at 298 °K. 
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average of about 3 cal mole" 1 with the accurate 
values [20]. Equation (28) gives H 2ai23) —H 2(298 ) = 
5837. From the enthalpy of the pertinent /? zirconium 
hydrides as formulated in section 6 (eq (92)), 
2[d(H° n23 -H° 298 )/dx] = 78l5. Using eq (27) and these 
values, eq (26) gives 



(H? +5 -#° 2 ) 298 =-41,478. 



(29) 



Equation (29) has an uncertainty of at least 1,000 cal 
mole -1 , but the additional significant figures serve 
to maintain consistency in many of the subsequently 
derived values. 

An equation equivalent to eq (26) in integral form 
may be derived which is more general because it 
indicates the energy difference between two zirco- 
nium hydrides where both x and T may differ finitely. 
If we consider a system of 1 g-atom of zirconium 
and 2 of hydrogen, and arbitrarily define the reference 
state as the two free elements in their standard 
states at a temperature of 298 °K, the increase in 
enthalpy when the temperature or chemical com- 
position or both change may be called the "intrinsic 
enthalpy" {H'). 1 The process is 

Zr(298°K)+H 2 (g)(298°K) 

= ZrH x (T) + (l-^)H 2 (g)(T), (30) 



and 



H'=H°+(l-lx)Hl-H° Ui 



"^2(298)- (31) 



From eq (13) 

1J-298 E 1 1(298) 2 • LL± 2 (298) 

= (m +s -WU+ix(Hl +s -H° 2 ) 2 



(32) 



Adding eqs (31) and (32), 

+hx{m +h -H° 2 ) 2fs8 +{l-hx){Hl-Hl (2m ). (33) 

This result may be applied to two compositions 
ZrHi and ZrH' x ' at T and T" respectively. Thus for 

ZrH I <(T')+Kx"-x')H 2 (g)(r)=ZrH I -,(T") (34) 

eq (33) leads to 

AH=(H°-H° 2 £ +i \n-(H°-H° 2 <r S) )*, 

+ W-x'){m^-H 2 ) 2 , & 

-W-x')(H° 2iTn -H° 2(29S) ). (35) 

In_ equation (26), (33), or (35), (H 2 -H 2i2 , S) ) and 
(H 2 +S —H 2 ) 2gs can be evaluated from eqs (28) and 
(29) respectively. Partial differentiation of eq (33) 
with respect to x and substitution from eq (6) show 
that 



{m'lbx) T =^{H 2 -H° 2 ). 



(36) 



With the help of eq (32), eq (5) may be put into 
convenient form for relating the total and partial 
molal enthalpy and free-energy functions in the 
form in which they are expressed in this paper, 
as well as immediately giving the standard heat and 
free energy of formation of ZrH^ at temperature T. 
In these two cases eq (5) becomes respectively 



and 



H =H x +lxH 2 

6°=6 1 +ix0 2 . 



(37) 
(38) 



Subtracting eq (37) from eq (32) (and remembering eq (13)) gives 



&HflT={H 1 -miT+\z(H % -JB®IT={H°-B& 

-hx{Hl-Hl {2m )IT+{Hr h -H\) 2 JT+^^^ (39) 

whereas subtracting eq (38) from eq (32) (and remembering eq (15)) gives 



AGf/T=(G-G° 1 )/T+hx(G 2 



-MG°2-H° 2(298) )/T+(Hr 8 - 



-m) 2 JT+\x(H« + >-Hl) 2 JT. (40) 



4. Alpha-Zirconium 

Pure metallic zirconium undergoes a first-order 
transformation from the hexagonal (a) to the high- 
temperature body-centered-cubic (0) lattice. The 
transformation temperature has been given values 
from 1,135 °K [21] to 1,143 °K [22, 23]. The value 
adopted in this paper is: 

Transformation temperature of Zr= 1,136 °K. (41) 



7 In an earlier paper [17] the author denned an equivalent quantity and called 
ti the "adjusted heat content." 



Measurements of the low-temperature heat ca- 
pacity of zirconium have been reported by Todd 
(53 to 297 °K) [24]; Skinner and Johnston (14 to 
298 °K) [25]; Estermann, Friedberg, and Goldman 
(1.8 to 4.2 °K) [26]; Wolcott (1.2 to 20 °K) [27]; and 
Burk, Estermann, and Friedberg (20 to 200 °K) 
[28]. The heat capacities of Todd and those of 
Burk, Estermann, and Friedberg are in close agree- 
ment, while those obtained above 130 °K by Skinner 
and Johnston are about 1 percent lower. The 
enthalpy at high temperatures has been measured 
by Mixter and Dana (273 to 373 °K) [29]; Jaeger and 
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Veenstra (294 to 1,074 °K) [30]; Coughlin and King 
(298 to 1,371 °K) [31]; Redmond and Lones (273 
to 1,309 °K) [32]; and Douglas and Victor (273 to 
1,173 °K) [16]. Furukawa and Reilly [33] have 
carefully analyzed the low-temperature data and have 
derived a table of thermodynamic functions from 
to 300 °K whose heat capacities approach those of 
Douglas and Victor at and above 300 °K. The 
heat-capacity curve of Douglas and Victor joins 
smoothly with that of Todd, but not with that of 
Skinner and Johnston. At 250 °K Furukawa and 
Keilly's tabJe gives (in the notation of sec. 2) 



and 



(^)i ( 250) = 5.969 
OS°) 1(250) =8.219. 



(42) 
(43) 



Using a sample of zirconium containing a total of 
0.09 percent by weight of impurities (0.015% Hf, 
0.03% Fe, 0.02% C, 0.005% N, and 0.006% O), 
Douglas and Victor [16] precisely measured the 
enthalpy, relative to 273 °K, at ioO-deg intervals 
up to 1,073 °K, and also at 1,123° and at 5-deg 
intervals from 1,143 to 1,173 °K. The following 
equation (in the notation of sec. 2) was derived to 
fit their data for the a form and to satisfy eq (42): 

H?-fl r H298.i5)=4.0504J f +5.3113 (lO" 3 ) T 2 

-4.4945 (10" 6 ) T 3 + 1.6785 (10~ 9 ) T 4 - 1573.9. (44) 

The mean observed values differ from this equation 
by the amounts shown in table 2. 

Table 2. Deviation of eq (44) from observed values of enthalpy 
of a- zirconium 





■^l 0- ^l°(273.2) 


Mean observed- 


Temperature 


calculated 




Mean observed 


enthalpy 


°K 


cal mole- 1 


cal mole- 1 


273. 2 


(0) 


(0) 


373. 2 


630 


+ 1 


473.2 


1284 


-9 


573. 2 


1989 


+8 


673. 2 


2686 


-1 


773. 2 


3409 





873.2 


4147 


_•> 


973.2 


4923 


+9 


1073. 2 


5708 


— 7 


1123. 2 


6133 


+ 1 



For the entropy of a-zirconium, eq (43) and (44) 
give 

Si=9.32645 log T+ 1.06226 (10" 2 ) T 
-6.74175(l(r 6 )r 2 +2.238 (10~ 9 ) T*-16A14:5. (45) 

The thermodynamic properties of ^-zirconium will 
be evaluated in section 7 in connection with the f$ field 
of the Zr-H system, of which it forms a part. 

5. The a and («+6) Phase Fields 

The thermodynamic properties of the a and 
(a +5) phase fields as derived in this paper are 
closely related, and will be discussed together. 
Other parts of the Zr-H system will be similarly 
treated in the following two sections. 



5.1. Choice of the a-Phase Eutectoid Point 

Values for the eutectoid temperature, where the 
a, f} and 8 phases of fixed compositions coexist in 
equilibrium, are in fair agreement as reported by 
various investigators. Vaughan and Bridge [10] 
gave 833°±10 °K, but the values most commonly 
reported have been closer to that of Ells and McQuil- 
lan [6], which is 820°±2 °K. In an enthalpy 
investigation [16], four samples transformed some- 
where between 823 and 833 °K, and one sample 
somewhere between 813 and 823 °K. The exact 
temperature is probably uncertain by a few degrees 
owing to some sluggishness in the (all solid) trans- 
formation, but for the purposes of the calculations 
of this paper the following value has been adopted: 



Eutectoid temperature = 820 °K. 



(46) 



Another constant evaluated empirically here is the 
eutectoid composition of the a phase. As indicated 
in figure 1, the a/(a+5) and a/(a-{-p) boundaries 
of different investigators [4,6,7,8,9] are in fairly good 
agreement. On the basis of these results the composition 
at 820°K selected here is 



z a (eutectoid) =0.0650. 



(47) 



The equilibrium properties of the (a +5) field are 
of course the appropriate combinations of the 
properties of the component phases along the two 
phase-field boundaries, which if known can be used 
to derive the two-phase-field properties. This pro- 
cedure has been followed for the (a+ff) field (sec. 7). 
If, however, more complete data are available in 
the two-phase field, these may be used to derive 
its properties even without knowing the boundaries 
of the field. In the present case some hydrogen- 
pressure data are available for the higher tempera- 
tures of the field, but not at the lower temperatures 
near room temperature, mainly because the equilib- 
rium pressures are negligibly small. 

5.2. Enthalpy and Related Properties 

The enthalpy measurements on five well-spaced 
compositions in the (a+8) field [16,17] show good 
precision and good consistency, and will be used as 
the basis for extending other properties below the 
eutectoid temperature. These data show not only 
the linearity with x of isotherms expected for a 
two-phase field [16,17], but show also other evidence, 
soon to be presented, that phase equilibrium was 
approximately maintained during the rather rapid 
cooling involved in the measurements. 8 repre- 
senting an isotherm of enthalpy by 



(H°-H° 29S . 15 )° +s =A+Bx, 



(48) 



the values of A and B (referred to the experimental 
base temperature of 273.15 °K) were determined 



8 Even if the eutectoid temperature is really as low as 820 °K (eq (46)), the 
(a+8) enthalpy isotherm at 823 °K is not to be interpreted as invalid, since the 
samples were not heated above 823 °K in these measurements and the /3 phase 
could easily have failed to appear. 
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from the data by the least-square method, and are 
given in table 3 referred to the experimental base 
temperature 273.15 °K — i.e., for the equation 



(H°-H° 273 . 15 )«+*=A'+B'x. 



(49) 



Table 3. Isotherms of enthalpy in the {a + 5) field of Zr-H 
(In cal mole- 1 of ZrH x ) 



T 


A' in eq (49) 


B' in eq (49) 


°K 






273. 15 








373. 15 


637.5 


114.0 


473. 15 


1305. 6 


356.5 


573. 15 


2041.8 


694.9 


673. 15 


2849. 4 


1039. 7 


773. 15 


3794. 2 


1410.0 


823. 15 


4277. 


1705. 8 



The following equations were derived to represent 
A and B in eq (48) as functions of temperature 

(T=273 to 823) 9 : 

A=33.91211T-5.051385(10- 2 )T 2 +4.76757(10- 5 )T 3 
-1.603496(10- 8 )r 4 +1.373028(10 6 )2 7 - 1 

-9.2289104(10 7 )T" 2 - 10324.4; (50) 

B = 3.05464r+5.5805(10- 4 )T 2 -7.098S95(10- 5 )T 3 
+2.608S68(10- 7 ) r 4 -3.393270(10- 10 ) T 5 
+ 1.520416(10- 13 )7 76 -447.8. (51) 

From eqs (50) and (51), the constants in eqs (48) 
and (49) are related by 

A'=A+158.1; (52) 

B'=B+19.0. (53) 

From eqs (1), (6), (8), and (48) 

(C p )i a+8 =dA/dT (54) 

and 

(O p ) 2 a+8 =2dB/dT. (55) 

From eqs (26), (28), (29), (48), and (51) 

( J ff 2 « +5 -fl r 2)/r=-40404-0.550727 7 +8.483(10- 4 )T 2 

-1.419779(10- 4 )r 3 +5.217736(10~ 7 )T 4 

-6.78654(10- 10 )T 5 +3.040832(10- 13 )T 6 

+ 12000I 7 - 1 . (56) 

The same properties were next formulated in the 
a field. Apparently no measurements of enthalpy 
are available except for zirconium metal (sec. 4). 
However, the a phases are all low in hydrogen 
(fig. 1), and the laws of dilute solutions should be 
applicable to good approximation. This involves 
the assumption that the enthalpy and heat-capacity 

9 Equation (51) agrees exactly with table 3. However, owing to a small com- 
putation error, eq (50) differs slightly from the table above 473 °K (by an aver- 
age of 7 cal mole -1 , which is within the precision of the data). 



isotherms are linear with x. It will presently be 
shown that at 273.15 °K the a field is extremely 
narrow; consequently, it is convenient to assume 
that at this temperature a-Zr is in the (a + 5) field. 
At 820 °K eq (44) then gives 

(fl^-flSft«0.-o=3753A (57) 

while eqs (47), (48), (50), (51), (52), and (53) give 

(HUr-H 2 ^J%. 0650=4359.8. (58) 

Thus, using eqs (52) and (53), 



#^ a) -#™=3594.9+9316.4x. 



(59) 



By assuming the hydrogen vibration frequency in 
ZrH 2 found by Flotow and Osborne [15], 1190 cm -1 , 
the following empirical equations were derived: 

T=298-820: (<^)?= -5.9446 +2.7364 (10- 2 )T 

-1.185(10" 5 )T 2 (60) 

T=820-1136: (Cy?=4.836+4.5(10~ 3 )r (61) 

Since these are independent of x, by eqs (9) and 
(8) (C p )i is identical with the heat capacity of 
a-Zr (from eq (44)). When the latter and eq (60) 
(or (61)) are substituted into eq (5) and then inte- 
grated to satisfy eq (59), we have: 

T=298-820: (H o <" ) -H o 3 f£. + 1 $ ) )/T=4.05Q4: 

+5.3113(10- 3 )r-4.4945(10- 6 )T 2 +1.6785(10- 9 )T 3 

-1578.5/r+[8243.0/r-2.9723+6.841(10- 3 )T 
-1.975(10- 6 )T 2 ]x; 
7=820-1136: (H° (a) -H° 2 ( 9 a 8 i 5 d) )/T 

=4.0504+5.3113(10- 3 )T-4.4945(10- 6 )T 2 

+ 1.6785(10- 9 )T 3 -1578.5/r 
+[6577.4/T+2.418+1.125(10- 3 )T]x. ( 63 ) 

The application of eqs (26), (28), and (29) to eqs 
(62) or (63) gives respectively: 

T= 298-820: (H? -H° 2 )/T= -23005/ T 

-12.6085+1.3414(10- 2 )T-3.950(10- 6 )7 72 ; (64) 

T=820-1136: {H a 2 ~H° 2 )IT 

^-26336/r-l. 828+1. 982(10- 3 )r. (65) 

5.3. Hydrogen Activity and the a/(o:+5) Boundary 

The thermodynamic functions (G^-G^/T, giving 
(by eq (10)) the hydrogen equilibrium pressures, 



(62) 
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will now be formulated l'<>r the a and (a-\-8) fields. 
As discussed in section :>, these will be based on eq 
(23) and the pertinent (hernial properties. The 
invariant eutectoid pressure derived in section 6 is 
the (a +8) value at that temperature: 



^2(820)- 



: 1.0442(10- ')• 



(66) 



It is reasonable to assume that in the dilute a phases 
the activity of atomic hydrogen is proportional to 
its concentration at any temperature, and hence a% 
is proportional to x 2 , as has been observed [4]. Then 
fromeqs (10) and (66) 



[(G^-G° 2 )/T] m = -18.2174, 
and, using eq (47), 



(67) 



[(G a 2 -Gl)IT) m =9.l~A7 log z-7.3536. (68) 

Substitution of eq (56) into eq (11) and integrating 
to satisfy eq (67) , 

(GS+'—Ql)/T= -40,403.48/7+1.26796 log T 
-8.4829(10- 4 )T+6000/r 2 +7.098842(10- 6 )T 2 
-1.739232(10- 7 )T 3 +1.696622(10- 10 )? 74 

-6.08162(10- l4 )T 5 +22.0491. (69) 

The similar use of eqs (11), (64), and (68) gives 

T=298-820: (0J-fl f S)/2 T =9.1517 log x 

-23005/7+29.02815 log T-\ .34113(10- 2 )T 

+ 1.97447(10- 6 )T 2 -54.2116, (70) 

while eqs (11), (65), and (68) give 

r=820— 1136: (G%-G 2 )/T=9.1517 log a— 26336/T 
+4.20976 log 2 7 -1.98179(10- 3 )2 T + 14.1221. (71) 

The a/(a+8) boundary at 820 °K has been re- 
quired to have a composition of # a =0.0650, but no 
other restrictions on it have been explicitly imposed. 
The boundary compositions at other temperatures 
can be found by equating either the hydrogen 
activities (eqs (69) and (70)) or the enthalpies (eqs 
(48), (50)/ (51), (62)) in the two fields. > The results 
by these two methods are compared in table 4. 10 
The agreement for the two methods above 400 °K 
is rather striking, and gives strong support to the 
assumption made earlier that in the a field a 2 is 
proportional to x 2 . If this assumption is correct, 
the agreement in the table indicates that phase 
equilibrium, especially with respect to the a/ (a +8) 
boundary, was approximately maintained in the 
enthalpy measurements in the (a + 8) field, a result 

i" The values In t he second column of the table are sensitive to the power of x 
assumed. With the assumption of the 820 °K value as before, the second column 

was recalculated with a a proportional to x 1 - 75 , but the disagreement with the 

third column was much worse. It may be noted, however, that the extent of 
disagreement is Independent ol the value assumed for the eutectoid hydrogen 
activity, for t lie same value was assumed in deriving eqs (69) and (70), and cancels 
when the two equations are set equal to each other. 



Table 4. 



The compositions of the al(a-\-h 
methods 



boundary by two 





x a at boundary 


T 








Equality of 


Equality of 




hydrogen 


enthalpy 




activity 




°K 






298. 15 


0.0000 15 


0. 0006 


300 


.000016 


. 0006 


350 


.00010 


.0010 


400 


. 00043 


. 0010 


450 


. 0013 


. 0013 


500 


.0031 


. 0025 


550 


.0064 


. 0053 


600 


.0114 


.0102 


650 


.DISS 


.0178 


700 


. 0288 


. 0282 


750 


.0418 


.0414 


800 


. 0578 


. 0576 


820 


.0650 


. 0650 



in accord with the observation of Schwartz and 
Mallett [9] that hydrogen diffuses in zirconium 
metal very rapidly at these temperatures. Aside 
from the small uncertainty in the value selected for 
820 °K, the values in the second column of table 4 
are probably more reliable than those in the third 
column; this is certainly true at the lowest tempera- 
tures, where the failure of the a phase to reach the 
final equilibrium compositions would have caused 
but little error in eq (56) and hence in eq (69). In 
this paper, however, the properties of the (a-\-8) 
field depending on the a/(a-\-8) boundary are calcu- 
lated from the third column of the table in order to 
avoid small discrepancies with other properties 
calculated from the thermal data. These discrep- 
ancies would have been minor, because the two 
columns of table 4 never differ by more than about 
0.001 in x. For the same reason, no effort was 
mad* 4 to adjust the (a + 5) enthalpy values to make 
them consistent with the second-column boundary 
compositions. 

5.4. Entropy 

The absolute entropies of the various hydride com- 
positions are calculated in this paper from those of 
zirconium metal (eq (45)). &£ can be calculated from 
eq (4) by substituting either eqs (64) and (70) or eqs 
(65) and (71), depending on the temperature. Spe- 
cifically, at 820 °K #2=38.285, and so we get 

Sj (M0) =13.3187-9.1516 log x. (72) 

Since from eq (45) 

S?(82o) = 16.1723, (73) 

eqs (72) and (73) may be substituted into eq (7) and 
the integration performed from x=0 to x ) giving 

S , ^? ) = 16.1723+8.6466a;— 4.5758a; log x. (74) 

The entropy may be formulated from the enthalpy 
by performing the integration (of eq (2)) to satisfy 
eq (74). When eq (63) is used for the enthalpy, the 
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result is 

r=820-1136: S 0(a) = (9.32645+5.5677x) log T 
+ [0.0106226+2.25(10- 3 )x]T-6.74175(10- 6 )7 72 
+ 2.238(10- 9 )7 73 -16.4145-9.4216;r-4.5758;rlogz. 

(75) 

Equation (74) gives (in view of eqs (46) and (47)) 

O0(a) O0(a + 5) —17 0874 (7(\} 

°820, 0.0650 ° 820, 0.0650 1 I .UCW±. \*V) 

At this temperature eqs (4), (56), and (69) and the 
value of ($2)320 give 



= 5.5207. 



(77) 



The substitution from eqs (76) and (77) into eq (5) 
then sives 



S? ( t&) = 16.9080. 



(78) 



(C P )i +s and (C p )% +S as functions of temperature were 
then evaluated from eqs (50, 54) and (51, 55) re- 
spectively, and when each of these functions was 
substituted into the equation for the partial molal 
properties analogous to eq (2) and the subsequent 
integration was required to satisfy eq (78) or (77) 
respectively, the partial molal and (from eq (5)) the 
total entropies in the (a+5) field were completely 
determined. The results for the partial molal 
entropies are 

S? +5 =78.08602 log r-0.1010277r+7.151355(10- 5 )T 2 
-2.13799467(10- 8 )2 73 +6.86514(10 5 )T- 2 
-6.15260693(10 7 )7 7 - 3 -164.9841; 



(79) 



^ +5 =14.06722 log r+2.2322(10- 3 )T 

-2.1296685(10- 4 )T 2 +6.95698132(10- 7 )7 73 
-8.483175(10- 10 )T 4 +3.6489984(10- 13 )2 75 

-29.4258. 



(80) 



5.5. Some Partial Molal Properties With Respect 
to Zirconium 

In the a field the partial molal relative free energy 
and relative enthalpy of zirconium may be found 
readily from those for hydrogen by integrating the 
Gibbs-Duhem equation (eq (9)), since zirconium in 
its standard state is a terminal composition (x = 0) 
of this field. Equations (70) and (71) both give the 
same result for the relative free energy, 



(fl f ?-0!)/r=-1.9872aj, 



(81) 



and eqs (64) and (65) give the same result for the 
relative enthalpy, 



(H?-H° 1 )/T=0. 



(82) » 



In the (a + 8) field the partial molal relative free 
energy of zirconium is simply the values given by 
eq (81) along the a/(a+8) boundary: 



(G? +8 -G < DIT=-l.9872x a . 



(83) 



The values of x a are given in table 4, where, as 
explained earlier, those in the third column were 
used to obtain the tabulated values of (G? +d —G° 1 )/T. 

The partial molal relative enthalpy of zirconium 
in the (a+8) field may be obtained from eq (4), 
which in this case is 

(H a 1 + *-H°)/T=(Gr 8 -G 1 )/T+Sr d -S 1 , (84) 

by substituting from eqs (45), (79), and (83). For 
the zirconium-hydrogen system this property may 
be derived more directly from the enthalpies alone 
by what is equivalent to the "method of intercepts." 
From eq (8) 



H^HV-xdH'/dx. 



(S5) 



By writing the same equation for 273 °K and the 
(a+8) field and subtracting from eq (85), we may 
write 

(H 1 -HV = (H«-H%r 8) )-x [^ (H°-m% +8) )~] T 

-(i??-i/; (273) ) + (^? +5 -^i)273. (86) 

It was stated earlier, however, that a-Zr will be 
assumed to lie in the (a + 8) field at 273 °K so that 

(Hr 8 -m) 27 *=o. (87) 

From eqs (48), (49), (52), and (53), 

(HT* +8) -m<r 8) )-x [^ (sss+w — ss& + «)] p =i58, 

(88) 
and from eq (44), 

£f; (298) -i?; (273) =154. (89) 

Subtracting eq (89) from the sum of eqs (86), (87), 
and (88), 



(H-H°)/T= 



H 



(H°-H^ +li ) 



dx 



I- 



(#»-i7? (298) )+4 



'}■ 



(90) 



For the (a + 8) field at T=298, eq (90) gives the value 
of eq (25). Equation (90) is applicable to any phase 
field and any composition lying inside the (a+d) 
field at 298 °K. In the special case of the (a+8) 
field, substitution from eq (48) gives the simple result: 



11 Equation (82) follows also from eqs (11) and (81). Both equations (81) and 
(82) are consistent with the assumption made earlier that the a phases are ideally 
dilute solutions of atomic hydrogen in a-zirconium. 
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(Hr 8 -m)/T={A-(m-H° U298} )+4]/T. (91) 

Values calculated from eqs (<S4) and (91) did not differ 
by more than 0.006, and were averaged for tabulation. 

6. The (8 (Hydrogen-Rich), (p+S), and 5 Phase 
Fields 

Having formulated thermodynamic properties for 
the (a + 5) field, those formulated for the 0, £+5, and 
6 fields should be consistent not only with the most 
reliable data for these fields but also with a single set 
of 0/(0+5), 5/(0+5), and 5/ (a + 5) boundaries. There- 
fore the properties of these three fields are considered 
together. 

An isotherm of a 2 at 1073.2 °K for the field has 
already been formulated in section 3 (eq (23)). As 
the hydrogen concentration (x) increases, the second 
term in this equation becomes relatively important 
(several percent of the first term) in the region x= 
0.5 to 0.6; other properties such as the heat capacity 
and the heat of hydriding (II2 — H2), at similar tem- 
peratures and the same general range of concentra- 
tions, show similar changes. This fact makes it 
convenient to treat the wide field separately in two 
parts. Somewhat arbitrarily, the field has been 
divided in this paper by the composition selected as 
that of the eutectoid phase, ^=0.57. The proper- 
ties of the phases richer in hydrogen will be con- 
sidered quantitatively in the present section, and 
those poorer in hydrogen, in section 7. The most 
appropriate dividing composition probably shifts 
with temperature, but the exact choice is probably 
not of major importance because all the important 
thermodynamic properties are made to have con- 
tinuity at this composition. 

6.1. The Enthalpy and Other Properties (x^0.57) 

A simple equation was selected to fit the enthalpy 
measurements [16] in the field for the higher values 
of x: 

x ^0.57: (H°w-H°^)/T= (9.664 + 2. lUx)T 

+ 1500z-2521.1. (92) 

The values calculated from this equation are coin- 
pared in table 5 with the mean observed enthalpies. 12 
(The latter were actually measured relative to 273.15 
°K, but have been made relative to 298.15 °K by 
subtracting the quantities indicated by eqs (52) 
and (53).) 

If any of the samples was actually inside the 
(0+5) field when measured, its enthalpy would be 
found to be lower than if at the same temperature it 
were in the field. Measurements were made on 
the compositions x= 0.556 and 0.701 also at 873.2 °K, 
but, contrary to figure 1, these enthalpies gave 
definite evidence (poor precision, large thermal 
hysteresis, and large negative deviations from eq 
(92)) that for some reason the field had not been 



Table o. Deviation from eq (92) of the mean observed 
enthalpies 

(For each temperature and composition the first value is the mean observed 
(/j°W— i/oj+fj) and the value in parentheses is the observed — calculated, 
both in cal/molc of ZrIT x .) 13 



T 


X- 


=0.556 


X- 


-0.701 


X = 0.990 


X= 1.071 


°K 

923.2 


8291 


±6(-44) 


8876 


±0(4-30) 






973.2 


8839 


(-39) 


9423 


(+24) 






1073.2 


9971 


(+8) 


10528 


(+13) 




[11906(-16)] 


1123.2 










[12211(— 27)] 


12520(0) 


1173.2 


11040 


(-9) 


11640 


(+8) 


12828(0) 


13118(0) 



reached. By similar considerations the observed 
values in table 5 for x=0.999 at T= 1123.2 and for 
x= 1.071 at T= 1123.2 and 1073.2 would be expected 
to fall below eq (92) . The fact that the last value is 
only slightly lower is evidence that the preceding 
value (for x= 1.071 and T= 1123.2) probably corre- 
sponds to the field, despite the 0/(0+5) boundary 
selected in figure 1. Most likely, small percentages 
of impurities in each sample, while corrected for 
additively, probably invalidated the observed en- 
thalpy appreciably only if the impurities shifted this 
phase-field boundary enough to place the point 
erroneously in a two-phase field. The doubtful 
values mentioned were not used in deriving eq (92). 
It may be noted from table 5 that the fit is poorer 
and with a trend for the two lowest temperatures 
tabulated (up to 0.5% discrepancy), but these four 
points were considered insufficient justification for 
modifying eq (92). 

With the help of eqs (28), (29), and (92), eq (26) 
gives: 

^0.57: (/7?-H5)/T=-36508.3/T-2.372 

-2.678(10- 4 ) T+ 12000/ T\ (93) 

When eq (93) is substituted into eq (11), and the 
subsequent integration is made to satisfy eq (24), 
there is obtained: 

£^0.570: (#|-#2)/r=6.2231 log x 

+4.5758 log (l + 1.8493/ ) - 97 ) 

+ 11.6100-36508.3/T+5.4617 log T 

+2.677(10" 4 )r+6000/T 2 . (94) 

Though eq (24) was derived from data up to about 
35=1, where the field terminates at 1073 °K, eq (94) 
is assumed to hold up to 1200 °K and also for the 
slightly greater values of x (up to about 1.1) which 
lie in the field above 1073 °K. An empirical equa- 
tion was derived for the standard entropy of hydrogen 
gas from eq (28) (using eqs (1) and (2)). Since eq 
(28) is not quite exact, the integrated equation was 
adjusted to give the best fit for T^820: 

T^820: S° 2 = 15.33532 log r+5.356(10- 4 )r 

-6000/T 2 -6.8288. (95) 



" The composition x =0.556 in table 5, while outside the range to which eq (92) 
has been restricted, is too close toj=0.57 to invalidate appreciably the differences 
"observed minus calculated" noted in the table. 



13 The tolerances attached to tw r o of the values in the table represent not pre- 
cision or estimated accuracy, but the uncertainty due to hysteresis. For these 
two points the temperature was approached both by heating and by cooling (in 
two separate series of runs) . The enthalpies found were higher in the former than 
in the latter series by twice the stated tolerance, and the mean of the two series 
has been tabulated. 
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Substitution from eqs (93), (94), and (95) into eq (4) 
then gives: 

z^0.57: "^ = -6.2231 logo; 

-4.5758 log (1 + 1.8493 2 6 - 97 ) + 9.8736 log T 

-20.8108. (96) 

Substituting from eqs (44) and (92) into eq (90), 

2^0.57: (0?--iZ?)/7 7 =--942.6/T+5.6136 

-5.3113(10" 3 ) T+4.4945(10- 6 ) T 2 - 1 .6785(10- 9 ) T\ 

(97) 

Other properties to be derived for the hydrogen-rich 
part of the (3 field will depend on how the (£+5) 
properties are formulated, because these, together 
with the /3 properties, will determine the unique 
composition of the ft eutectoid phase with which they 
are consistent. 

6.2. Procedure Giving Thermodynamic Consistency 

Of the twelve thermodynamic functions derived 
and tabulated in this paper (tables 8-14), five must 
be continuous at each phase-field boundary: H, S, 
G, 6ri, and G 2 . Both relative-enthalpy and hydrogen- 
activity data have been obtained in the (0+5) field. 
The enthalpy data were obtained by Douglas and 
Victor [16] for the compositions #=0.701, 0.999, and 
1.071, but the values show serious inconsistencies 
between two of the sample compositions, 14 and 
therefore were given no weight in deriving properties 
in this field. Hydrogen-activity values in the 
(jS+5) field were reported as follows: by Edwards 
et al. [5] for x= 1.222 at 973, 1,023, 1,073, 1,098, 
1,123, 1,133, 1,138, and 1,148 °K; by LaGrange et al. 
[8] for several compositions at 843, 889, 919, 967, 
1,022, and 1,073 °K. The values of the two sets 
of workers, which disagree somewhat at the lower 
temperatures, show excellent agreement as well as 
smooth temperature dependence at and above 1,023 
°K, and were used to define the adopted values in 
this higher temperature range. With the complete 
formulation of the hydrogen activity (or G 2 ^ 8 ) from 
the eutectoid temperature (820 °K) to 1,200 °K, 
soon to be given, the /3/(/3+<5) boundary was so 
determined as to make G% and Gi +d equal along the 
boundary, and the values of the other four of the 
above five properties in the (0+5) field were deter- 
mined so that along the boundary they would be 
equal at each temperature to the corresponding 
/3-field properties. 

Representing the boundary composition at tem- 
perature Thy xp, the relative enthalpy in the (0+5) 
field at x and T can be found by integrating eq (26) 
between x and x&\ 



The complete determination of (H%+ 8 —HI) at the 
lower temperatures down to 820 °K will then com- 
plete the determination of the (0+5) properties 
and the 0/(0+6) boundary. In particular, the 
composition of the eutectoid phase will lie on this 
boundary, and the composition of the eutectoid 
8 phase will be that resulting from simultaneous 
solution of eqs (48) and (98) at 820 °K, the eutectoid 
temperature. It was decided to formulate thermo- 
dynamic properties of the Zr-H system only up 
to #=1.25, partly because there are fewer over- 
lapping data for compositions richer in hydrogen 
and those that do exist show some major incon- 
sistencies. If the eutectoid 8 phase has a value of 
x less than 1.25, then as the temperature rises the 
composition #=1.25 passes from the (a + 5) field into 
the 5 field at some temperature (TO, and from the 
8 into the (0+5) field at a higher temperature (T 2 ). 
For thermodynamic consistency, it must be possible 
to formulate the 5-phase properties so that the five 
properties mentioned above agree with those of the 
(a + <5) and (0+5) fields at the eutectoid 5 composition 
and temperature (820 °K), with the (a + 5) properties 
at #=1.25 and T u and with the (0+5) properties at 
2 = 1.25 and T 2 . 

The (0+6) hydrogen-activity data above_1073 °K 
mentioned above indicate a value of {Hl +b —H° 2 ) 
essentially independent of temperature (approxi- 
mately — 51600), and in a preliminary trial the same 
value was assumed down to the eutectoid temper- 
ature. This assumption led, in the manner outlined 
above, to eutectoid compositions of approximately 
££ = 0.54 and 2 5 =0.99. The next step taken was to 
formulate an enthalpy function for the 6 phases 
whose temperature derivative corresponded to a 
hydrogen vibration frequency of between 1050 and 
1200 cm" 1 [15] and which at the same time satisfied 
the (a+6) enthalpy equation (eq (48)) at the eutec- 
toid 6 point and also on the 6/ (a + 5) boundary at 
673 °K. With 2=1.25, this equation was then 
solved simultaneously with eq (48) to determine T h 
and with eq (98) to determine T 2 . To test the con- 
tinuity of the other four properties at the points T x 
and T 2 , entropy is an insensitive criterion, but G 2 
proves entirely adequate. The 6-phase enthalpy 
function mentioned above and eq (26) gave a value 
of approximately (Hl—W 2 ) = — 42900, and on sub- 
stitutino; this value, eq (11) integrated to give for 
m-G° 2 )/T] T2th25 -[(GI-G 2 )/T] T ,. 2 , about 0.7 



(H°« + »-H° 2 £.lP)IT-- 



less than [(Gi +8 -G^)/T] T ^. 25 -[(G^-GI)JTl Tl>1 ^ 
(The discrepancy of 0.7 was approximately 15% of 
either of the two differences.) 

This serious discrepancy was practically eliminated 
by making two changes. The first change was to 
preserve the fit to the hydrogen-activity data at the 
higher temperatures but to assume slightly higher 
values of (W +8 —Hl) at the lowest temperatures 

= (ff°w-tf?&+^ [(2-2,)/2T]. (98) 



1 4 The discrepancies are much larger than those found in several tests of heating- 
and-cooling hysteresis, and seem to be due to different locations of the 0/08+5) 
boundary caused by different impurities in the two samples. 
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( n cur 820 °K). The final equations are 

(Hi+*-m)iT= 

[-51678.2 +3.1277(10- 10 ) (1200- T) 5 ]/T (99) 

and (from eq (11) after substituting from eq (99)) 

(G^ 8 -Gl)/T = 726592.5/T+ 7466.81126 log T 

-5.40465827+2.2519409(10-*) T 2 

-6.2553915(10- 7 )T 3 +7.81924(10- n )r 4 

-19434.0399. (100) 

The additional eutectoid compositions (at 820 °K) 
were then found in the ways outlined above. Simul- 
taneous solution of eqs (94) and (100) gave 



^(eutectoid) =0.570. 



(101) 



Using this value, simultaneous solution of eqs (48) 
and (98) (after substitution from eqs (50) and (51) 
into the former, and from eqs (28), (29), (92), and 
(99) into the latter) gave 



ara (eutectoid) = 1.171. 



(102) 



The effect of increasing the value of j 5 (eutectoid) 
from 0.99 to 1.171 was to bring the boundary tem- 
peratures T x and T 2 (at #=1.25) much closer to- 
gether with the result that the discrepancy in the 
interval of (G 2 —Gl) x= 1.25 referred to above was now 
greatly reduced, but still had the same sign, and a 
magnitude of 0.14. To repeat the above change on 
an increased scab 4 and thus to increase the eutectoid 
composition x$ still further would have increased 
the discrepancies of the p/(a + fi) boundary (sec. 7) 
but apparently ^without eliminating the present 
discrepancy in (G\— Gl)i .25. This discrepancy was, 
however, practically eliminated by a second change 
which, by lowering (JJl — H 2 ), makes a small increase 
in the temperature coefficient of (G 2 —G 2 )/ r r'u) the 
<5 field. This was accomplished by replacing the 
5-phase enthalpy function by the simple equation 

z=1.17 to 1.4, 7=790 to 850: (H 5 -H£+ 5 15 )/T 

= 13.115- (4692.5+ 19.0s) IT. (103) 

With the use of eqs (103), (28), and (29), eq (26) 
now gives 

z=1.17 to 1.4, 7 = 790 to 850: (H* 2 -Ht)/T 
= -39546.3/T-6.660 

-2.678(10- 4 )r+ 12000/ T\ (104) 

Equation (104) gives approximately (H 8 2 —H 2 °) = 
— 45200, which incidentally is much closer than the 
preliminary value —42900 to that derived by 
Gulbransen and Andrew [4] from their data, —45750. 
The 0/(0+5) boundary on which eqs (94) and (100) 
give the same values of G 2 is that shown in figure 1, 
and for convenience may be represented (to within 
±0.001 in x) by the empirical equation 
7 = 820 to 1200 °K: x = - 1. 3995 + 3. 5169(10" 3 )7 7 

-1.2(10- 6 )7 72 +[1.19853(10- 15 ) 



Equation (101), when substituted into eq (94) (or 
(100)), gives for the eutectoid activity of hydrogen 
the value of eq (66) given earlier, 1.0442(10" 4 ) atm. 15 
The heat capacity which eq (103) gives corre- 
sponds (by eqs (1) and (6)) to 



(cyi=o, 



(106) 



and this deserves some comment. Though eq (103) 
was derived to fit only a narrow range of temperatures 
and compositions, no claim can be made that eq (106) 
is highly accurate and hence that the addition of 
hydrogen to the 8 phases produces an intrinsic addi- 
tion to the heat capacity by the vibration of the 
added hydrogen which is exactly offset by a reduction 
of the intrinsic contribution by the vibration of the 
zirconium. It seems not unreasonable, however, 
that a considerable compensation of this kind may 
actually exist. Equation (80) gives negative values 
for S" at the lowest temperatures such as 298 °K, 
which means that the addition of hydrogen to a- 
zirconium to form the coexistent 5-phases lowers 
the entropy at these temperatures, and hence pre- 
sumably the heat capacit ics at all lower temperatures. 
In these cases it is reasonable to assume that a major 
result of adding the hydrogen has been to stiffen the 
zirconium lattice, and that this relative effect is 
dependent on composition as well as on temperature. 

6.3. Other Properties (./ ^0.57) 

With the value of eq (101), further properties 
may be derived in the f$ field. From eqs (5), (77), 
and (78), 

S° 8 ^ ) 57o=18.4814. (107) 

Equations (48), (50), and (51) give 

(H°-H° 298 . 15 ) So. 570=5040.5, (108) 

and eq (92) gives 

(H°W-ffSa+i«) 820,0.570=7260.5, (109) 

whence 

A£ 820 =2.7073 for ZrH . 570 (a + 5) -> ZrH o . 5 7o(0). (HO) 

Adding eqs (107) and (110), 

S°a&o.57o=21.1887. (HI) 

For the composition #=0.570 eqs (92) and (111) 
readily give the entropy: 

S£ ( $ o =25.06629 log T-51 .8497. (112) 

After substituting from eqs (96) and (112), eq (7) 
can be used to derive the /3 entropies for #>0.570. 

15 The values indicated by the pressure measurements of Gulbransen and 
Andrew 14] and of Ells and McQuillan [6] seem to fall within a range of ±10 
percent of this value. 

-2.87092(10- 18 )r+1.6655(10- 21 )r 2 ] [1200-T] 6 . (105) 
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This involves the integral 

-2.2879 | log (1 + 1.8493 z 697 ) dx. Because the 

J 0. 570 

expansion of the integral about x=0 does not con- 
verge rapidly for values of x approaching unity, the 
following empirical equation was derived: 

z=0.57 to l.l:-2.2879 | log (1 + 1.8493 x 697 ) dx 

J 0. 570 

= 0.4202-2. 12581:z+3.72544x 2 

-2.36774^+0. 18558^. (113) 

The general entropy equation then found was 

£^0.57:£ 0( ^ = (22.25230+4.93681;z) log T— 46.7017 
-11.17991x+3.72544x 2 -2.36774x 3 
+0.18558x 4 -3. 111552: log x. (114) 

S^(x^0.57) was calculated from eqs (5), (96), and 
(114). 

Since from eqs (47) and (101) ZrH . 65o(«) an d 
ZrH . 57 o(/5) are coexistent phases at 820 °K, eq (81) 
(with x a =0. 0650) evaluates Gi for both phases: 



[(^-^)/r] 820 , .57o=-0.1292. 



(115) 



When eq (97) is substituted into eq (11) and the 
integration is made to satisfy eq (115), there is 
obtained for the one composition x= 0.570 

[G?-0?)/T]o. 570= -942.6/I 7 - 12.92588 log T 
+5.3113(10- 3 )2 7 -2.24725(10- 6 )7 72 +5.595(10- 10 )T 3 

+35.5312. (116) 



Equations (94) and (116) may be substituted into 
the Gibbs-Duhem equation (eq (9)), the integration 
from £=0.570 to higher x giving (with the help of 
eq (113) in integration by parts) 

z^0.57:(G?-G'?)/T=35.9017 + 0.77449:z-3.72544£ 2 

+2.36774x 3 -0.18558x 4 -2.2879£ log (l + 1.8493x 6 - 97 ) 
-942.6/T-12.92588 log T+5.3113(10- 3 )T 

-2.24725(10- 6 )T 2 +5.595(10- 10 )T 3 . (117) 

6.4. Other (0+8) Properties 

The remaining (13+5) properties may be readily 
derived from foregoing relationships. Si +8 may be 
found from eqs (4), (95), (99), and (100). The 
entropy was found from that at the 0/(0+5) bound- 
ary (S° x { ^, calculated from eqs (105) and (114)) 
and eq (7) in the form 



Sow+«=S»(«+JSg + «(x-^) 



(118) 



S?+ 5 was calculated by substituting S 0((i+6) and Sl +8 
into eq (5). (G$ +s —G[)/T, which is independent of 
x, is the value given by eq (117) on the boundary 
and at the temperature in question. (Hj +S —H%)/T 
was calculated from eq (4) from Sl +8 , (G^-G^/T, 
and S?(eq (45)). Cl (p+i) is the total derivative of eq 
(98) with respect to temperature (with x but not 
Xj3 independent of temperature); by eq (5) this may 
be broken up into the partial molal properties: 



(C p )i +i =[b(H°^-IP^ 1 p)/bT] x ,+WIP' 



m i2g8) )]/dT 



-M(^r 5 -^)-(^ + 5 -ff2°)298+(^-^( 29 8))][^^]; 



(119) 



(c P )^=d[(H^-m)+(m-m i2g8) )]/dT. (120) 

(The quantities needed in eqs (119) and (120) can 
be obtained from eqs (28), (29), (92), (99), and (105).) 

6.5. Other 8 Properties 

The remaining properties in the small part of the 
8 field considered here (from 7"= 790 to 850, and 
from x= 1.171 to 1.25) were derived as follows. 
From eqs (90), (103), and (44), 

3=1.17 to 1.4, T=790 to 850: (H'-H^/T 

= 9.0646-5.3113(10- 3 )T+4.4945(10- 6 )T 2 

-1.6785(10- 9 )T 3 -3114.0/T. (121) 

The following boundary points were selected as 
giving the best overall continuity of properties at 
the points: 



x=1.25:8/(a+8) boundary, 
Ti=791. 5; 8/(0+5) boundary, T 2 =845.6. (122) 
At these points eqs (69) and (100) give respectively 

[(^ +5 -^)/r] 79L5L25= ,_20.0576; 

[(#I +5 -<?2)/TW 6il . 25 = -16.3872. (123) 

However, the integration between these two tem- 
peratures of eq (11), after substitution froin eq (104) 
gives a corresponding increment of (G$ — G$)/T 



which is 0.0203 Jess than that of eq (123); but by 
adjusting the integration constant so that the dis- 
crepancy has at the two points the same magnitude 
but opposite signs, the same equation then gives 

[(#S-G^)/T] 8 2o,i.25= -18.0684. (124) 

From eqs (67) and (102) 

[(#£-<?g)/T] 8 2<u.i7i=-18.2174. (125) 
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By making the reasonable assumption that in this 
small composition range a* is proportional to some 
constant power of 2, eqs (124) and (12r>) <rive: 

x=1.17to L25: [(Gl-Gl)/T] B20 

=5.2281 log x- 18.5751. (126) 

After substitution from eq (104), the general inte- 
gration of eq (11) so as to satisfy eq (126) then gives: 

x = 1.17 to 1.25, T=790 to 850: (G 8 2 -G%)/T 
= 5.2281 log Z-15.2606-39546.3/T+15.33532 log T 

+2.678(10" 4 )r+6000/T 2 . (127) 
From eqs (90), (103), and (44), 

3=1.17 to 1.4, T=790 to 850: (H*-H$/T 

=9.0646-5.3113(10- 3 )Z 7 +4.4945(l()- ,i )T 2 

-1.6785(10- 9 )T 3 — 3114.0/T. (128) 
From eqs (102) and (115) we have 



[(<?!-0!)/TW i.m= -0.1292. 



(129) 



Equations (126) and (129) may be substituted into 
the Gibbs-Duhem equation (eq (9)), the integration 
from 2=1.171 to higher x giving 

2=1.17 to 1.25: [(Gl-G1)/T] 820 =l. 1997-1. 13526 x. 

(130) 

When eq (128) is substituted into eq (11) and the 
latter is then integrated so as to satisfy eq (130), 
there is obtained 

2=1.17 to 1.25, T=790 to 850: (Gi-G^)/T 

-1.13526 2+62.6621-31 14.0/T- 20.87215 log T 

+ 5.3113(10- 3 )T-2.24725(10- 6 )T 2 + 5.595(10- 10 )2 73 . 

(131) 

The partial molal entropies can be derived from 
eq (4)^ Substitution from eqs (45), (128), and (131) 
gives Si, and substitution from eqs (95), (104), and 
(127) gives S 8 2 : 

2=1.17 to 1.25, T=790 to 850: Si 

=30.19860 log T+ 1.13526 2-70.0120; (132) 



2=1.17 to 1.25, T=790 to 850: S* 2 

= 1.7718-5.2281 log 2. 



(133) 



The substitution of eqs (132) and (133) into eq (5) 
gives the entropy: 



2=1.17 to 1.25, T=790 to S50: S o(d) 

= 30.19860 log T+2.02116 x 
-2.61405 2 log 2-70.0120. 



(134) 



For the composition 2=1.25 the discrepancy at 
each adopted phase-field boundary between the final 
equations for the <5 and the two-phase fields is given 
in table 6 for each of the five properties which should 
have identical values. (In each such case, the 
average value is tabulated.) Since continuity was 
required at the boundary at 2=1.171, presumably 
the discrepancies would be small at the boundaries 
for all intervening values of x. 

Table 6. Discrepancies of formulated properties of ZrHi.25 
at the adopted phase-field boundaries 



Property 


"h" 
at 


-"(a+fi)" 
791.5 °K 


'>S"-"(fi+5)» 
at 845.6 K 


-(< 


s° 

o(or + 5)X 
298.15 )' 




+0 .005 
+ .002 
-.002 


+0 .007 
+ .008 
+ .001 




{Gx-G^jT 




- .006 


+ .006 




(G 2 -G 2 )/T 




+ .010 


- .010 



6.6. Discussion of the 8 Phase-Field Boundaries 

The 8/(a + 8) and 5/ (0+5) boundaries have been 
drawn in figure 1 through the adopted 8 eutectoid 
point (eq (102)) and the adopted temperatures for 
2=1.25 (eq (122)). Though no effort is made in 
this paper to derive consistent thermodynamic 
properties for higher values of 2, these two boundaries 
have been extended in figure 1 in a reasonable but 
empirical way. In particular, the ensemble of 
experimental points shown for the 5/(0+5) boundary 
shows wide discrepancies, 17 but the boundary drawn 
has been made to pass through the "hydrogen- 
pressure" values of Edwards et al. [5] and of LaGrange 
el al. [8]. The points of LaGrange et al. indicate 
some curvature in the boundary, but those of 
Edwards et al. (at seven temperatures from 973 to 
1,148 °K) show a vertical boundary within the 
experimental error. In a previous paper [17] the 
author criticized a vertical boundary on thermody- 
namic grounds which would lead one to predict that 
it has a positive finite slope. The same argument 
still holds, but may be reexamined somewhat more 
quantitatively. 

The hydrogen activity in a one-phase field b is of 
course a function of both temperature and compo- 
sition, and its variation with temperature along 
the boundary x b (T) separating the b and (6 + c) 
fields may be written as a total derivative: 

d In a 2 /(lT=(d In a 2 /dT) x +(d In a,/dx) T /(dT/(Ix b ), 

(135) 
or in the present case 

d In a b 2 +c /dT=(d In al/dT) Xb + (d In a f 2 ) /dx) T /(dT/dx b ). 

(136) 



'6 For the eutectoid 8 phase Or— 1.171, 7*=820), eq (134) should of course give 
the same value of S° as eqs (5), (79), and (80) for the («+ &) field. The two values 
calculated were found to differ by only 0.002. However, this is a check of thermo- 
dynamic consistency only. 

" Relatively accentuated for higher values of x in figure 1 compared with the 
more usual plots of the phase diagram as temperature versus atomic percentage 
of hydrogen. 
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Substituting for the first two derivatives from eqs 
(10) and (11), eq (136) becomes 



(din al/dx) T 
= (RT*)-\dT/dx b )[m-m)-(H$ +c -H° 2 )]. 



(137) 



If the thermodynamic properties as formulated are 
consistent, regardless of how near to correct they 
are, they must satisfy eq (137) along every phase- 
field boundary. For example, the equation may be 
applied to the 8/(a+8) and 8/(0+8) boundaries at 
#=1.25. In the former case the bracketed quantity 
in eq (137) is -3006 (from eqs (56) and (104)), and 
in the latter case it is +4577 (from eqs (99) and 
(104)). Since, as expected, (d In a 2 8 /(>x) T is positive 
in both cases (eq (127)), the slopes of the 8/(a+8) 
and 8/(0+8) boundaries at this composition will 
thus be negative and positive respectively, as they 
are in figure 1. 

The measured 5-phase isotherms of hydrogen 
pressure [3,5,8] are in genera] considerably steeper 
than indicated by eq (127), according to which a\ is 
at constant temperature approximately proportional 
to x 114 (up to x= 1.25) . It may be noted furthermore 
that the isotherms cited usually show a more gradual 
change in slope, 18 as the hydrogen content changes 
and the sample crosses a phase-field boundary, 
when the latter is contiguous with the 8 field than 
when it is contiguous with some other one-phase 
field. If, as seems to be the case, the heats of 
hydriding do not change rapidly with temperature, 
eq (137) indicates that at values ol x considerably 
greater than 1.25 the far steeper 8/(0+8) boundary 
indicates a much greater value of (d In a 8 2 /dx) T than 
at lower values of x, as experiment indicates. By 
the same reasoning, the 8/(a+8) boundary would be 
expected to become much steeper as x increases. 
This boundary in figure 1 was drawn to pass through 
the points indicated by thermal-expansion studies 
from 293 to 673 °K [12,14], and does become 
somewhat steeper. If this boundary were assumed 
to be an equilibrium one, it would be possible to 
draw some conclusions about the 5-field hydrogen 
activities at temperatures too low for direct accurate 
measurements. 

A sample such as ZrHi 2 o shows the interesting be- 
havior of being in the one-phase 8 field at some 
temperatures, but in a two-phase field at higher or 
lower temperatures. According to equations derived 
earlier (5, 50, 51, 54, 55, 103, 119, and 120), the 
equilibrium heat capacity increases discontinuously 
as the sample passes from the 5-field into either 
contiguous two-phase field, because a new phase of 
different energy starts to form at a finite rate. A 
simple proof may be given that the heat capacity is 
necessarily greater, near the boundary, in the two- 
phase than in the one-phase field. If ' and " indi- 
cate respectively the one- and two-phase properties, 
at the boundary G" = G', S" = S', and so also 
dG"/dT=dG'/dT since S=-dG/dT. But since in- 



side the two-phase field G"<G' ', then the G" (T) 
curve for the sample must be more negatively curved 
than the partly metastable G' (T) curve, or d 2 G"/dT 2 
<d 2 G f /dT 2 . Hence dS"/dT>dS'/dT, and since dS/ 
dT=C v /T, it follows that C 9 ">C 9 '. 

7. The j3 (Hydrogen-Poor) and (a+0) Phase 
Fields 

Unlike the thermodynamic properties in the (a + 8) 
and (0+8) fields derived in sections 5 and 6, those 
in the (a+0) field were derived from the properties 
of the two adjacent one-phase fields using phase-field 
boundaries thermodynamically consistent with these 
properties. This consistency was imposed by find- 
ing two boundaries such that at each temperature 

GS = G$ m fm& G? = G% (eqs (18) and (22)). The a-field 
properties were formulated in section 5, and the 
0-field properties for x^ 0.570, in section 6. How- 
ever, since it is assumed in this paper that x^ 0.570 
in all of the (a+0) field (eq (101) and fig. 1), the 
pertinent /3-field properties are those for the lower 
range of x ) and are derived in the present section. 
This includes /3-zirconium (x=0) from its trans- 
formation temperature, 1,136°, up to 1,200 °K. 

7.1. Discussion of Thermodynamic Consistency 

The thermodynamic data in the presently con- 
sidered ranges of temperature and composition which 
were given greatest weight are as follows. Figure 1 
shows, in addition to the a/ (a+0) and 0j(a + 0) 
boundaries finally arrived at in this paper, various 
points on these boundaries indicated by three ex- 
perimental investigations [6, 8, 10]. The points 
from the two investigations of hydrogen pressure 
[6, 8] showed general agreement, but the X-ray 
measurements [10], which did not involve the a/ 
(a + 0) boundary, indicated a much higher curve for 
the 0/(a+0) boundary. The compositions of the 
eutectoid phase reported by these investigators are 
equivalent to the following values of x: Ells and 
McQuillan [6], 0.47; LaGrange et al. [8], 0.50; 
Vaughan and Bridge [10], 0.72. Douglas and 
Victor's [16] measurements of the enthalpy of 
Zr(jS) and of ZrH 324 (0), when reduced to a basis of 
298.15 °K by eqs (52) and (53), give the following 
mean values: 



(TT°(P) —TI°( a + V\ —7014. 

1^1173.2 ^298.15 ^1 = *0±*±, 



ytt 1173. 2 



-HM)o. 3 24=9701; (138) 



18 This curvature of the 5-field isotherms increases the difficulty of locating the 
compositions of the breaks accurately, and suggests that some of the values of 
x selected in the references cited may be somewhat too high. 



in addition, their mean observed values of (i?° (/3) 
-H ^) o.32 4 were 8836 at T= 1073.2 and 8004 ± 2 19 
at T=973.2. They actually measured the enthalpy 
of Zt(0) at five temperatures over the very short 
temperature range 1,153 to 1,173 °K, the values 
varying linearly with temperature with a heat capa- 
city averaging 9.0. Skinner's linear heat capacity- 
temperature equation, which is based on his enthalpy 
measurements of Zv(0) from the transformation 



19 See footnote 13 (table 5) . 
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temperature to 1,N00 °K [22], gives 6.73 in the same 
temperature range and extrapolates to 6.46 at 1,000 
°K. 

Development of 0-phase properties for £^0.570 
began with the assumption of the eutectoid tem- 
perature a- and /3-compositions already adopted 
(eqs (47) and (101)), the adopted transformation 
temperature of zirconium (eq (41)), and the basic 
relations in sections 2 and 3. In the first trial an 
enthalpy equation for the region was formulated 
in the following way. It was required to satisfy 
the values of eq (138), to agree at £=0.57 in H and 
(dll/dx) r with the equation already derived for 
#^0.57 (eq (92)), and to correspond to a heat 
capacity independent of temperature but equaling 
Skinner's 1,000 °K value for x=0 (6.46) and Douglas 
and Victor's mean value for x= 0.324 (8.48) . 20 
In the manner outlined later, points on the a/(a-\~l3) 
and 0/(a-f-j8) boundaries at several temperatures 
were then found by imposing the requirement that 
the two partial mola] free energies be equal. 

The /S/(a-(-jS) boundary so found lay considerably 
above that indicated by the hydrogen -pressure 
studies of Ells and McQuillan [6] and LaGrange 
et al. [8] and considerably below that of Vaughn n 
and Bridge [10] (except at x=0, where there is 
general agreement). While the X-ray technique 
is less sensitive than the hydrogen-pressure method, 
it seems reasonable to conclude that unless Vaughan 
and Bridge made a gross misinterpretation of their 
data their samples actually did correspond to a 
boundary in the neighborhood of what they re- 
ported and that some factor such as unknown 
impurities accounts for the large differences from 
the reported hydrogen-pressure curves. 21 Neverthe- 
less the latter two curves, although they show some 
differences at the lower temperatures, are in sub- 
stantial agreement, and seem much more credible. 
In the author's derivation of the two boundaries as 
outlined above, the single requirement that the 
hydrogen activities be equal (G%=Q§) gives of 
course an infinite set of curves, one of which approxi- 
mated the experimentally observed "hydrogen pres- 
sure" boundaries, but at intermediate temperatures 
with zirconium activities of approximately 0.94 on 
the a boundary and 0.97 on the fi boundary. Just 
how sensitive these zirconium activities are to shifts 
under the experimental conditions of the hydrogen- 
pressure measurements is not apparent. 

It is appropriate of course to examine critically the 
assumptions and basic data used by the author in 
deriving the boundaries. The a-field activities of 
hydrogen and zirconium given by eqs (70), (71), and 
(81), which are related by the exact Gibbs-Duhem 
equation, 22 are not only reasonable but are supported 
by the agreement in table 4 below 820 °K. An 
independent entropy calculation showing excellent 



-' Douglas and Victor's value of 9.0 for 0-Zr, which is greater than they found 
for ZrHo.334, is not precise because of the short temperature range, and may be 
too high If the small amount of oxygen and nitrogen (0.065 atom %) raised the 
transformation temperature above 1,153 °K [61. 

2i Ells and McQuillan [f>] concluded, e.g., that oxygen as an impurity raises 
this boundary. 

22 The boundaries arc Insensitive to the eutectoid hydrogen pressure used as an 
integration constant . 



agreement will be made later in this section, and 
supports the basic hydrogen-activity isotherm as- 
sumed (eq (23)). The boundaries are insensitive at 
the lower temperatures to the heat capacity assumed 
for jS-zirconium. The enthalpy measurements on 
ZrH 324 were of high precision (±5 cal mole -1 ), the 
sample was of fairly high purity, and there are good 
reasons for believing that in tin 4 temperature range 
cited above (973 to 1,173 °K) the sample was entirely 
in the one-phase field and hence not subject to the 
uncertainties caused by phase reproportionation. 
Nor did it seem justified to modify eq (92) for the 
"hydrogen-rich" /3 enthalpies. In fact, it seems clear 
that enthalpies corrected additively for small amounts 
of impurities may still be seriously in error in multi- 
phase fields owing to improper phase compositions, 
but reliable in one-phase fields. In several attempts 
it was found impossible to reformulate the 0-phase 
enthalpies for x ^0.570 with the restrictions men- 
tioned above and still obtain a credible approximately 
monotonic function which had an appreciable effect 
on the derived a/(a + (3) and fi/(a+0) boundaries. 
Furthermore, the choice of a somewhat lower value 
of the eutectoid ft composition would do some violence 
to the fit which led to eq (101) without leading to a 
much better fit to the experimental "hydrogen- 
pressure" p/(a + j3) boundaries. This is because the 
latter suggest a boundary slope near the eutectoid 
temperature which is several times less steep than 
that required by eq (137). 

7.2. The (3 Properties (a^O.57) 

In order to derive a thermodynamically consistent 
set of properties in the range of temperature and 
composition under consideration (IT>82(), x<0.570), 
it was decided to give approximately equal weight 
to the hydrogen-pressure results [6, 8] and the en- 
thalpy data [16] by making such changes in the latter 
as would lead to a derived (3/(a + /3) boundary ap- 
proximately halfway between the experimental 
' ' hydrogen-pressure" boundary and that derived 
above. A new £-field enthalpy function was derived 
by the procedure described earlier except that the 
heat capacity of Zr(/3) was raised from (5.46 to 7.76 
and that of ZrH 0324 (/3) was raised from 8.48 to 9.00: 23 

x^ 0.570: (H°^-Hl^ + J ) )/T=7.7^]+2AUx 

-f-90.047x 4 -126.382r 5 
+ (4990.1x-109834.5x 4 + 147541. Or 5 - 1 793. S)/ J. (139) 

From the temperature derivative of eq (139), eqs 
(6) and (8) give the partial molal heat capacities: 

zSO.570: (<7 p )?=7.763-270.141z 4 +505.528x 5 ; (140) 

zg0.570: (("gg=4.2SS+720.376^-1263.N2.r l . (141) 

From eqs (44), (90), and (139): 

23 By retaining the values of eq (138), the enthalpies of ZrHo.324 below 1,173 °K 
were varied by many times what seems the probable uncertainty of the observed 

values. However, this procedure had the desired effeel of increasing (_dH^ldx)T at 
the lower temperatures near the phase-field boundary. 
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x^ 0.570: (m-H\)IT 

= (-215.3+329503. 5a; 4 -590166.4a; 5 )/T 
+3.7126-270.141a: 4 +505.528a; 5 -5.3113(10- 3 )T 

+4.4945(10- 6 )T 2 -1.6785(10- 9 )T 3 . (142) 

From eqs (26), (28), (29), and (139): 

a;;g 0.570: (^-iT 2 )/r=-2.372+720.376a; 3 -1263.82a; 4 
+ (-29528.1-878676a; 3 + 1475416a; 4 )/r 

-2.678(10- 4 )T+12000/T 2 . (143) 

By substituting from eq (143) in eq (11). integrating 
to satisfy eq (24), and adding a small correction 
term "A" (a function of x only) which will be pres- 
ently evaluated (eq (151)), there is obtained 

x^ 0.570: (Ol-GDlT 

-6.22313 log x+4.5758 log (1 + 1.8493a; 6 - 97 ) 
+ 5. 1053+ 5845.8474a; 3 - 10194.2788a; 4 
+ (-29528.1 -878676x 3 +1475416x 4 )/T 
+ (5.46177- 1658.73778a; 3 +2910.07193a; 4 ) log T 
+2.678(10- 4 )r+6000/T 2 +A. (144) 

On substitution from eqs (95), (143), and (144), eq 
(4) gives (with log (1 + 1.8493a; 6 - 97 ) expanded to one 
term) : 

»£ 0.570: Sg=-6.22313 log x -3. 675a; 6 - 97 

-14.3061-5125.4714a; 3 +8930.4588a: 4 

+ (9.87355 + 1658.73778a; 3 -2910.07193a; 4 ) log T-A. 

(145) 

The entropy will next be found. From eq (45) 

the entropy of a-zirconium at the transformation 
temperature is 

S!fil,.o=18.7293. (146) 

Equations (44) and (139) give for the heat of trans- 
formation of zirconium 



(H o ^-H°^) imtZ=0 =937.1. 



(147) 



Adding the corresponding entropy of transformation 
to eq (146), 

#?i£<Uo= 19.5542. (148) 

The entropy at 1,136 °K and any value of x up to 
0.570, found from eq (7) by integrating from a;=0 
to x after substituting from eqs (145) and (148), is 

x^ 0.570: S S=19.5542+9.28232a;-3.11156a;loga; 

-0.23a; 7 - 97 -7.17464a; 4 +3.9083a; 5 -i J Adx. (149) 

When, using the heat capacity from eq (139), eq (2) 
is integrated to satisfy eq (149), the result is 



x^ 0.570: S O( < 3) = (17.87508+4.93677a; 

+207.3422z 4 -291. 0072x 5 ) log T-35.0610 
-5.80139a;-3.1 1156a; log x-0.23x 797 
-640.6838a; 4 +893.0459a; 5 -^ | Adx. (150) 



The way in which eq (150) was derived indicates no 
value of A other than zero. With A=0, the equation 
gives for x =0.5 70 values of entropy which are higher 
than those of eq (112) by 0.016 at all temperatures. 
To make the two equations agree, A was arbitrarily 
chosen to be 



A= -0.13 log (z/0.570) 



(161) 



in eqs (144), (145), (149), and (150). Agreement 
with the value of eq (148) is maintained, and for 
£=0.57 eqs (94) and (144) still agree. Since eq (151) 
changes the calculated values of hydrogen activity 
by only 3 percent lor £=0.2 and by less for higher 
values of x, the agreement with the observed values 
as tabulated in table 1 is still good and within the 
experimental error. 

Equations (112) and (150) (with A=0) both 
depend on the entropy of zirconium metal at 820 °K, 
but were otherwise derived by different paths of 
composition- and temperature-change from this 
common point. The small disagreement of 0.016 
noted above is thus a good check on some of the 
data entering the derivations, particularly on the 
basic hydrogen-activity equation (eq (23)) since 
this equation enters into the derivation of eq (150) 
only. However, some of the thermal data assumed 
affect both equations to comparable extents. Spe- 
cifically, the close agreement is no real check on the 
P enthalpy functions (eqs (139) and (143)), nor on 
the exact choice of the eutectoid j3 composition (eq 
(101)), for to these it is rather insensitive. 

When eq (144) (with A from eq (151)) is substi- 
tuted into the Gibbs-Duhem relation (eq (9)) and 
the integration is required to satisfy eq (116), there 
is obtained 

x ^0.570: (#?-~<7 1 )/T=22.3595--1.32311z 
-2192.19277a; 4 +4077.7115a; 5 - 1.607a; 7 - 97 + 1.6a; 14 - 94 

+ (-215. 3 + 329503. 5x 4 -590166.4x- 5 )/r 

+ (-8.54863 + 622.02667a- 4 - 1164.02877a; 5 ) log T 

+ 5.3113(10- 3 )T-2.24725(10- 6 )T 2 

+ 5.595(10- 10 )r 3 . (152) 

For a;=0 and T=1136, this equation gives +0.004. 
This value should be zero, since 1,136 °K has been 
assumed to be the transformation temperature of 
zirconium. S\ for x ^0.570 can readily be found 
from eq (4) or (5) using equations already given (45, 
142, and 152). 
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7.3. Derived (a + /3) Boundaries 

Points x a and X3 on the a/(a + 0) and j8/(a+£) 
boundaries as determined by the formulations for the 
a and fields given in this paper were then calcu- 
lated at several temperatures, and are listed in table 
7. This was done by requiring the respective partial 
molal free energies to be equal at the boundary com- 
positions at the temperature in question (eqs (71), 
(81), (144), (151), and (152)). The compositions at 
the eutectoid temperature had already been deter- 
mined (eqs (47) and (101)). The boundaries meet 
at the transformation temperature of zirconium. 
The solid a/(a-\-ff) and £/(«+£) boundaries in figure 
I were drawn through the points in table 7, except 
that they were allowed to deviate slightly from the 
1,000 °K points, meeting as straight lines at T= 1136 
and x=0. LaGrange et al.f8] have pointed out that 
the two boundaries are required thermodynamically 
to meet at this point with a difference in slope equal 
to the heat of transformation of zirconium divided 
by RT 2 . In the notation of the present paper this 
relation is 



/dxg dx 8 \ 
\dT dTj im , x 



\ irn 



/11 



(153) 



and is true provided the activity of diatomic hydro- 
gen is proportional to x 2 in both solid solutions. 
Despite eq (23), this is a reasonable assumption for 
very dilute solutions, and the two boundaries have 
been drawn in figure 1 so as to meet with slopes 
satisfying eq (153) with the heat of transformation 
of eq (147). At the higher temperatures the pro- 
cedure of this section in defining the boundaries be- 
comes increasingly sensitive to errors in the data, and 
no calculations of this nature were attempted above 
1,000 °K. 



Table 



'. — Calculated boundary compositions of flu 
(a -f #) phase field 



T 


Xa 


Xf} 


°K 






820 


0. 0650 


0.570 


850 


. 0696 


.535 


865. 3 


.0717 


.500 


873.2 


.0727 


.472 


900 


.0747 


.355 


950 


.0745 


.259 


9m. 6 


.0741 


.250 


1000 


.0708 


.202 


1136 


.0000 


.000 



7.4. The (a+/3) Properties 

The values of (S? + ^-G?)/T and (G% +(i r G 2 )/T are 
of course the values at the boundaries given by the 
equations just cited. The enthalpy function and the 



entropy, (H o(a+ 0-fl8£& l) )/T and S 0(a+ ^\ are given 
by eq (16) after substituting from eqs (63), (75), 
(139), (150), and (151) and table 7. S% + ? and S? +fl 
were calculated from eqs (6) and (5) respectively. 
(H? +6 -m)IT and (m + ^-m\/T were then obtained 
from eq (4). The latter function was also calculated 
from eq (26), which takes the form 



(n% +(i -H°)iT= 

[2/T(x 3 -xM(H^-H^ +8) )r,-(Il 0{a) - 

+(m^-m) 2 jT-(m-m i29H) )/T. (ua) 



ZT0(a+6)\ 

1 ' 2'.)S )xa 



The values calculated by the two methods differed 
(owing to small computational inconsistencies) by 
0.007 on the average, and the mean was tabulated. 24 
The heat capacity in the (a+P) field was formu- 
lated by substituting the a and enthalpies (given by 
eqs (63) and (139) into eq (16) and then taking the 
total derivative with respect to temperature. For 
simplicity replacing (H ^-H° 2 £ + ») Xa by "H«" and 
{H om —WAl +b) ) X(i by "fl»" and omitting the sub- 
scripts x a and sty, the result is 



(JO (a + 



w= /aH^ 



dT A 



~\_Xff-x a J L (>T + \ dx " x s -x a ) (IT J 

The phase-field boundary slopes in eq (155) may be 
replaced in terms of more common thermodynamic 
functions. Since 



G%+»=G" 2 -- 



we have 



■-GS, 



(156) 



djGS+'-GD/T b(G2-GZ)lT 



dT 



dT 



b(Gj-GI)ITdx b 



dx 



dT' 



(157) 



where b is either a or 0. dxJdT and dxp/dT may be 
substituted into eq (155) from eqs (157). If the 
temperature derivatives appearing in eq (157) are 
replaced from eq (11), some simplification is achieved 
by substituting for the resulting partial molal 
enthalpies from eqs (26) and (154). Equation (155) 
finally becomes 



La*-*J L ST "•" V dx Xf>-x a Jr <>x\ T /J 



-mm<^- H i^yK(^)] <*> 



The partial molal heat capacities can be found from 
eq (158) by using eq (5). 



24 As a further check, the values of (f/? + *-Hf)/T and (i^-J^)/Tagreed 
approximately with those estimated from eq (11). 
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The procedure outlined in this section is generally 
applicable to determining the phase-field boundaries 
and the thermodynamic properties of two-phase 
fields from the properties of the adjoining one-phase 
fields. In the present case, there is considerable 
uncertainty in the (a+/3) properties because of 
uncertainty in the exact locations and shapes of the 
boundaries, and this is particularly true of the heat 



capacity, which eq (155) shows to be highly de- 
pendent on the slopes of the boundaries. However, 
as noted earlier, small changes in the one-phase 
properties often have large effects on the calculated 
boundaries. For this reason, due consideration of 
direct evidence as to the positions of the boundaries 
can serve in a treatment of the present type to define 
the properties of the neighboring one-phase fields 
more accurately. 



Table 8. Thermodynamic functions for zirconium, Zr 
T in deg K, thermodynamic functions in cal (deg K)-i mole -i. Subscript 1 refers to Zr (a) . Subscript 2 refers to PT 2 (g) . 
(C P )i=C p o 8i = S° ^ 2 =-(G 2 -(?2)/r=oo 























T 


Phases 
present 


H°-H°m.\5 


c° 

p 


S° 


G°-H 298.15 


H°-H i 


Gi-Ol 


lh-Hl 


(C,)i 




T 


T 


T 


T 


T 




298.15 


a 





6.197 


9.290 


-9.290 








-86.12 


1.16 


300 


a 


.038 


6 .205 


9.329 


-9.291 








-85.62 


1.20 


350 


a 


.934 


6.404 


10 .301 


-9.367 








-74.13 


2.18 


400 


a 


1.628 


6.572 


11.167 


-9.539 








-65.39 


3.10 


450 


a 


2.186 


6.712 


11.949 


-9.763 








-58.49 


3.97 


500 


a 


2.644 


6.830 


12 .663 


-10.019 








-52.90 


4.78 


550 


a 


3.030 


6.931 


13.319 


-10.289 








-48.25 


5.52 


600 


a 


3 .359 


7 .020 


13 .926 


-10.567 








-44.32 


6.21 


650 


a 


3 .643 


7.102 


14 .491 


-10.848 








-40.95 


6.84 


709 


CY 


3.893 


7.182 


15 .020 


-11.127 








-38.02 


7.40 


750 


a 


4.115 


7.265 


15.519 


-11.404 








-35.44 


7.91 


791.5 


a 


4 .2S2 


7.340 


15.912 


-11.629 








-33.53 


8.29 


800 


a 


4.315 


7.357 


15 .990 


-11.675 








-33.16 


8.36 


820 


a 


4 .390 


7 .396 


16.172 


-11.782 








-32.32 


8.53 


845.6 


a 


4 .482 


7.451 


16.400 


-11.919 








-31.30 


8.64 


850 


a 


4.497 


7.461 


16 .439 


-11.942 








-31.13 


8.66 


865.3 


CL 


4.550 


7.496 


16 .573 


-12.023 








-30.55 


8.73 


891.5 


a 


4.637 


7.561 


16.797 


-12.160 








-29.60 


8.85 


900 


a. 


4 .665 


7 .584 


16 .869 


-12.204 








-29.31 


8.89 


950 


a 


4.822 


7.729 


17 .283 


-12.461 








-27.67 


9.11 


956.6 


a 


4.842 


7.751 


17 .337 


-12.494 








-27.46 


9.14 


1000 


a 


4.972 


7 .904 


17.684 


-12.712 








-26.18 


9.34 


1050 


a 


5.116 


8.111 


18.074 


-12.958 








-24.83 


9.56 


1080 .9 


a 


5.204 


8.258 


18.312 


-13.108 








-24.05 


9.70 


1100 


a 


5.258 


8.357 


18.457 


-13.199 








-23.59 


9.79 


1136 


a 


5.359 


8.560 


18.729 


-13.370 








-22.76 


9.95 


1136 





6.18 


7.76 


19.55 


-13.37 


0.82 





-28.66 


4.29 


1150 


8 


6.20 


7.76 


19.65 


-13.45 


.80 


-0.01 


-28.35 


4.29 


1200 





6.27 


7.76 


19.98 


-13.71 


.73 


-.04 


-27.29 


4.29 



Table 9. Thermodynamic functions for ZrH .25 
T in deg K, thermodynamic functions in cal (deg K)-i mole-i. Subscript 1 refers to Zr(«). Subscript 2 refers to H~2(g). 



T 


Phases 
present 


11 "298.15 


o% 


S° 


@ ^298.15 
f 


Hx—H° 


(C p )x 


Si 


Cn—Go 


H 2 -H° 


(C,)i 


Si 


G 2 -G% 




T 






T 


T 


T 


T 


298.15 


a+8 





6.54 


9.27 


-9.27 


0.02 


6.32 


9.31 


-0 .00 


-139.12 


1.71 


-0.29 


-107.62 


300 


a+8 


0.04 


6.54 


9.31 


-9.27 


.02 


6.32 


9.34 


-.00 


-138.29 


1.74 


-.28 


-106.76 


350 


a+8 


.98 


6.78 


10 .33 


-9.35 


.02 


6.42 


10.33 


-.00 


-119.20 


2.86 


+ .07 


-86.95 


400 


a+8 


1.73 


7.11 


11.26 


-9.53 


.02 


6.58 


11.19 


-.00 


-104.73 


4.24 


.54 


-72.02 


450 


a+8 


2.35 


7.51 


12.12 


-9.77 


.03 


6.81 


11.98 


-.00 


-93.32 


5.57 


1.12 


-60.36 


500 


a+8 


2.89 


7.94 


12.93 


-10.05 


.04 


7.12 


12.71 


-.00 


-84.07 


6.56 


1.76 


-51 .03 


550 


a+8 


3.37 


8.39 


13.71 


-10.35 


.08 


7.51 


13.41 


-.01 


-76.44 


7.05 


2.41 


-43.38 


600 


a+8 


3.80 


8.85 


14.46 


-10.66 


.13 


7.97 


14.08 


-.02 


-70.06 


7.03 


3.03 


-37.01 


650 


a+8 


4.21 


9.32 


15.19 


-10.98 


.21 


8.48 


14.74 


-.04 


-64.69 


6.75 


3.58 


-31 .62 


700 


a+8 


4.59 


9.86 


15.90 


-11.30 


.31 


9.01 


15.39 


-.06 


-60.09 


6.76 


4 .08 


-27.00 


750 


a+8 


4.97 


10.53 


16.60 


-11.63 


.43 


9.53 


16.03 


-.08 


-56.07 


8.00 


4.58 


-23.00 


800 


a+8 


5.34 


11.48 


17.31 


-11.97 


.55 


10.00 


16.66 


-.11 


-52.41 


11.86 


5.20 


-19.49 


820 


a+8 


5.50 


11.99 


17.60 


-12.10 


.61 


10.17 


16 .91 


-.13 


-50.98 


14.54 


5.52 


-18.22 


820 


a+8 


6.49 


9.99 


18.59 


-12.10 


.26 


7 .56 


16.56 


-.13 


-40.26 


19.43 


16.24 


-18.22 


850 


«+0 


6.64 


11.60 


18 .97 


-12.34 


.24 


6.88 


16.82 


-.14 


-38.15 


37.74 


17.19 


-16.81 


900 


ct+8 


7.18 


17.32 


19.90 


-12.73 


+ .07 


3.96 


17 .09 


-.15 


-31.22 


106.89 


22.53 


-14.80 


950 


a+8 


7.61 


14.82 


20.74 


-13.13 


-.08 


5.26 


17.35 


-.15 


-25.53 


76.49 


27.08 


-13.27 


956.6 


a+8 


7.66 


14.94 


20 .84 


-13.18 


-.10 


5.29 


17.39 


- .15 


-24.86 


77.14 


27.62 


-13.10 


956.6 


8 


7.66 


8.53 


20.84 


-13.18 


+ 1.23 


7.20 


18.72 


-.15 


-35.49 


10.61 


16.99 


-13.10 


1000 


8 


7.70 


8.53 


21.22 


-13.52 


1.15 


7.20 


19.04 


- .20 


-33.80 


10.61 


17 .46 


-11.56 


1050 





7.74 


8.53 


21 .64 


-13.90 


1.06 


7.20 


19.39 


-.25 


-32.03 


10.61 


17.98 


-9.96 


1100 


8 


7.77 


8.53 


22.03 


-14.26 


0.96 


7.20 


19.72 


- .30 


-30.42 


10.61 


18.47 


-8.50 


1136 





7.80 


8.53 


22.31 


-14.51 


.89 


7.20 


19.95 


-.33 


-29.35 


10.61 


18.82 


-7.54 


1150 





7.80 


8.53 


22.41 


-14.61 


.86 


7.20 


20 .04 


-.34 


-28.96 


10.61 


18.95 


-7.18 


1200 





7.84 


8.53 


22.77 


-14.94 


.76 


7.20 


20 .35 


-.38 


-27.61 


10.61 


19.40 


-5.98 
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8. Tables of Thermodynamic Functions 

The common integral and differential thermo- 
dynamic properties were calculated as described 
previously in this paper, and are listed in tables 
8-14 for zirconium and several compositions of 
zirconium hydride (ZrH x ). In addition, table 15 
gives the properties of ideal "normal" hydrogen gas 
(H 2 , 25% para and 75% ortho) [20], which are 
closely related to those of the Zr-H system. The 
tables cover evenly spaced compositions of ZrH* 
with x=0, 0.25, 0.50, 0.75, 1.00, and 1.25; in addi- 
tion, the composition £=0.57 is also included in the 
tabulation, parti}' because this composition is of 



special interest as the assumed /3-phase eutectoid 
composition, and partly because the /3-field properties 
are separately formulated on both sides of this 
composition. In conformity with the conventions 
used throughout this paper, temperatures are in 
degrees Kelvin and the thermodynamic properties 
are given in defined t hermochemica] calories per 
degree per mole of Zrll x (or H 2 , in table 15). The 
temperatures tabulated, which begin with 298.15 
°K, run every 50° from 300 to 1,200 °K except that 
the eutectoid temperature, 820 °K, and the trans- 
formation temperature, 1,130 °K, are included in 
all the tables because of some interest in these par- 
ticular isotherms. In addition, each table includes 



Table 10. Thermodynamic functions for ZrHo.so 
T in deg K, thermodynamic functions in cal (deg K)-' mole- 1 . Subscript 1 refers to Zr(o:). Subscript 2 refers to 112(g). 































T 


Phases 

present 


H°-HI, SM 


C° 


s° 


" '■'2A8.15 


Hx-Hf 


(Cp)i 


Si 


Gi-0f 


Ih-H- 


(C P h 


Si 


Gt-G$ 




T 






T 


T 






T 


T 






T 


298. 15 


«+5 





6. 75 


9. 23 


-9. 23 


0.02 


6. 32 


9. 31 


-0. 00 


-139. 12 


1.71 


-0.29 


-107.62 


300 


ck-H 


0.04 


6.76 


9. 28 


-9.23 


.02 


6. 32 


9.34 


-.00 


-13S.29 


1.74 


-.28 


-106.76 


350 


a+5 


1.03 


7.14 


10. 34 


-9.32 


. 02 


(1. 42 


10. 33 


-.00 


-119.20 


2. 86 


+ .07 


-so. 95 


400 


tt-N 


1.82 


7.64 


11.33 


-9.51 


. 02 


6. 58 


11. 19 


-.00 


-104.73 


4.24 


.54 


-72.02 


450 


a+5 


2.50 


8.21 


12.20 


-9.7(3 


.03 


6.81 


1 1 . 98 


-.00 


-93.32 


5.57 


1.12 


-60. 36 


500 


«+5 


3.10 


8.76 


13. 15 


-10. Of, 


.04 


7.12 


12. 71 


-.00 


-S4.07 


(1. 5(1 


1.76 


51.03 


550 


a+5 


3.64 


9. 27 


14.01 


-10. 3S 


.OS 


7.51 


13.41 


-.01 


-76. 44 


7. 05 


2.41 


-43. 38 


GOO 


a+5 


4.12 


9.73 


14.84 


-10.72 


.13 


7.97 


14. OS 


-.02 


-70.06 


7. 03 


3.03 


-37. 01 


650 


a+5 


4.57 


10. 17 


15. 64 


-11.06 


.21 


8.48 


14.74 


-.04 


-64. 69 


6. 75 


3.58 


-31.62 


700 


a+5 


4. 99 


10.70 


16.41 


-11.42 


.31 


9.01 


15. 39 


-.06 


-00.09 


6. 70 


4.08 


-27. 00 


750 


a+5 


5.40 


11.53 


17.17 


-11.78 


.43 


9. 53 


10.03 


-.08 


-56.07 


s. 00 


4.58 


-23. 00 


800 


a+5 


5. 82 


12. 97 


17.9(1 


-12.14 


.55 


10.00 


16.66 


-.11 


-52.41 


11.86 


5.20 


-19.49 


820 


a+5 


6.00 


13. 80 


18.29 


-12.28 


.61 


10.17 


16. 91 


-.13 


-50.98 


14. 54 


5. 52 


-18.22 


820 


a-r-0 


8.34 


12. 42 


20. 62 


-12.28 


. 2(1 


7.56 


L6.56 


-.13 


-40. 26 


19.43 


10.24 


-18.22 


850 


«+0 


8.54 


16.32 


21.12 


-12.59 


.24 


6.88 


L6.82 


-.14 


-38. 15 


37.74 


17.19 


-1(1. SI 


865.3 


a+0 


8. 72 


22. 57 


21.46 


-12.74 


. 22 


5.74 


16.94 


-.14 


-36. 72 


07. 33 


18.09 


-16.14 


865. 3 





8.72 


10.51 


21.41'. 


-12.74 


2.55 


(1. (is 


19.2(1 


-.14 


-40.02 


15.35 


8. 78 


-1(1.14 


900 





S. 79 


10. 51 


21.87 


13.0.x 


2.41 


(i.(iS 


19.52 


-.24 


-43.93 


15.35 


9.39 


-14.37 


950 





s. ss 


10.51 


22. 44 


-13.56 


2. 24 


(1. OS 


10. SS 


-.37 


-41.19 


15. 35 


10. 22 


-12.07 


1000 






8.96 


10.51 


22. 98 


-14.02 


2. 07 


6.68 


20. 23 


-.48 


-38. 72 


15.35 


11.00 


-10.02 


1050 


9.04 


10.51 


23. 49 


-14.46 


1.91 


6. 68 


20.55 


-.57 


-30.49 


15.35 


11.75 


-8. 19 


1100 


/3 


9. 10 


10.51 


23. 98 


-14.88 


1.75 


6. 68 


20.SC, 


—.66 


-34.4(1 


15.35 


12. 47 


-6.54 


1136 





9.15 


10. 51 


24. 32 


-15.17 


1.64 


6.68 


21.08 


-.71 


-33. 11 


15.35 


12.97 


-5.45 


1150 





9. in 


10.51 


24. 45 


-15.29 


L.59 


6.68 


21. 10 


-.73 


-32. 61 


15. 35 


13.16 


-5.05 


1200 





9. 22 


10.51 


24.90 


-15.68 


1.44 


6.68 


21. 14 


-.80 


-30.92 


15.35 


13. SI 


-3.70 



Table 11. Thermodynamic functions for ZrHo.57 
T in deg K, thermodynamic functions in cal (deg K)-i inole-*. Subscript 1 refers to Zr(«). Subscript 2 refers to 112(g). 































T 


Phases 

present 


"°-^2°9S.15 


C° 


S° 


G°-H° 

298.15 


Hi-Hf 


(C p )i 


Si 


Oi-Of 


lh-U% 


(C P h 


£2 


7h-<i° 2 




T 






T 


T 






T 


T 






T 


29S. 15 


a+5 





6.81 


9.22 


-9.22 


0.02 


6.32 


9.31 


-0.00 


-139.12 


1.71 


-0.29 


-107.62 


30(1 


a+5 


0.04 


6.82 


9.26 


-9. 22 


.02 


6.32 


9.34 


-.00 


-138.29 


1.74 


-.28 


-100.70 


350 


a+5 


1.04 


7.24 


10.34 


-9.31 


.02 


6.42 


10.33 


-.00 


-119.20 


2.86 


+.07 


-86.95 


400 


a+5 


1.85 


7.79 


11.35 


-9.50 


.02 


6.58 


11.19 


-.00 


-104.73 


4.24 


.54 


-72.02 


450 


a+5 


2.54 


8.40 


12.30 


-9.76 


.03 


6.81 


11.98 


-.00 


-93.32 


5.57 


1.12 


-60.36 


500 


a+5 


3.16 


8.99 


13.22 


-10.06 


.04 


7.12 


12.71 


-.00 


-84. 07 


6.56 


1.76 


—51. 03 


550 


a+5 


3.71 


9.52 


14.10 


-10.39 


.08 


7.51 


13.41 


-.01 


-76.44 


7.05 


2.41 


-43. 38 


000 


a+5 


4.22 


9.98 


14.95 


-10.73 


.13 


7.97 


14.08 


-.02 


-70.06 


7.03 


3.03 


-37.01 


650 


a+5 


4.68 


10. 40 


15.76 


-11.09 


.21 


8.48 


14.74 


-.04 


-64. 69 


6. 75 


3. 58 


—31. 62 


700 


a+5 


5.10 


10.94 


16.55 


-11.45 


.31 


9.01 


15. 39 


-.06 


-60.09 


6.76 


4.08 


-27.00 


750 


a+5 


5. 52 


11.81 


17.33 


-11.82 


.43 


9.53 


16.03 


-.08 


-56. 07 


8. 00 


4.58 


-23.00 


800 


a+5 


5.96 


13.38 


18.14 


-12. 18 


.55 


10.00 


16.66 


-.11 


-52.41 


11. so 


5.20 


-19.49 


820 


a+5 


6.15 


14.31 


18.48 


-12.34 


.61 


10.17 


16.91 


-.13 


-50.98 


14. 54 


5. 52 


- 18. 22 


820 





8.85 


10.89 


21.19 


-12.34 


2.21 


9.66 


18.51 


-.13 


-47. 10 


4.29 


9.41 


-18.22 


850 





8.93 


10.89 


21.58 


-12.65 


2.21 


9.66 


18.86 


-.21 


-45.53 


4. 29 


9.5(1 


-16.55 


900 





9.04 


10.89 


22.20 


-13.17 


2.20 


9.66 


19.41 


-.34 


-43. 16 


4. 29 


9.80 


-14.02 


950 





9.13 


10.89 


22.79 


-13.66 


2.19 


9. 66 


19.93 


-.45 


-41.04 


4.29 


10.03 


-11.74 


1000 





9. 22 


10.89 


23.35 


-14.13 


2.18 


9.66 


20.42 


-.56 


-39. 14 


4. 29 


10.2.-, 


-9.69 


1050 





9.30 


10. 89 


23.88 


-14.58 


2.15 


9.66 


20. 90 


-.67 


-37.41 


4. 29 


10. 46 


-7.82 


1100 





9.37 


10.89 


24.39 


-15.02 


2.12 


9.66 


21.35 


-.77 


-35.85 


4. 29 


LO.66 


-0. 12 


1136 





9.42 


10.89 


24.74 


-15.32 


2.09 


9.66 


21. 66 


-.84 


-34.80 


4. 29 


10.80 


-4.98 


1150 





9.44 


10.89 


24.87 


-15.43 


2.08 


9.66 


21.78 


-.86 


-34. 42 


4.29 


10. so 


-4.56 


1200 





9.50 


10.89 


25.33 


-15.84 


2.03 


9.66 


22. 19 


-.95 


-33. 11 


4. 29 


11.04 


-3. 12 
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the temperature at which the derived properties 
indicate that the particular composition crosses a 
phase-field boundary, with properties being tabu- 
lated for both the fields which meet at that boundary. 
Except for the properties of hydrogen gas and for 
those of a-zirconium which were not derived from 
those of the zirconium hydrides, the values of the 
properties have been rounded off to two decimal 
places as being more in keeping with their uncer- 
tainties in general. 
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Table 12. Thermodynamic functions for ZrHo.75 
T in (leg K, thermodynamic functions in cal (deg K)-i mole - *. Subscript 1 refers to Zr(«). Subscript 2 refers to H2(g). 































T 


Phases 
present 


" "298.18 


C° p 


S° 


G°-H° 

298.15 


H x -H° x 


(C p )i 


Si 


Gi-G° 


H 2 -H° 


(C p ) 2 


#2 


Ot-ai 




T 






T 


T 






T 


T 






T 


298.15 


a+8 





6.97 


9.20 


-9.20 


0.02 


6.32 


9.31 


-0.00 


-139.12 


1.71 


-0.29 


-107.62 


300 


a+8 


0.04 


6.98 


9.24 


-9.20 


.02 


6.32 


9.34 


-.00 


-138.29 


1.74 


-.28 


-106.76 


350 


a+8 


1.07 


7.49 


10.35 


-9.28 


.02 


6.42 


10 .33 


-.00 


-119.20 


2.86 


+ .07 


-86.95 


400 


a+8 


1.91 


8.17 


11.40 


-9.48 


.02 


6.58 


11.19 


-.00 


-104.73 


4.24 


.54 


-72.02 


450 


a+8 


2.65 


8.90 


12.40 


-9.75 


.03 


6.81 


11.98 


-.00 


-93.32 


5.57 


1.12 


-60.36 


500 


a+8 


3.31 


9.58 


13.37 


-10.06 


.04 


7.12 


12.71 


-.00 


-84.07 


6.56 


1.76 


-51 .03 


550 


a+8 


3.91 


10.16 


14.32 


-10.41 


.08 


7.51 


13.41 


-.01 


-76.44 


7.05 


2.41 


-43.38 


600 


a+8 


4.45 


10.61 


15.22 


-10.77 


.13 


7.97 


14 .08 


-.02 


-70.06 


7.03 


3.03 


-37.01 


650 


a+8 


4.94 


11.01 


16.08 


-11.15 


.21 


8.48 


14.74 


-.04 


-64.69 


6.75 


3.58 


-31.62 


700 


ct+8 


5.39 


11.55 


16.92 


-11.53 


.31 


9.01 


15.39 


-.06 


-60.09 


6.76 


4.08 


-27.00 


750 


a+8 


5.83 


12.53 


17.74 


-11.92 


.43 


9.53 


16.03 


-.08 


-56.07 


8.00 


4.58 


-23.00 


800 


a+8 


6.30 


14 .45 


18.61 


-12.31 


.55 


10.00 


16 .66 


-.11 


-52.41 


11.86 


5.20 


-19.49 


820 


a+8 


6.51 


15.62 


18.98 


-12.47 


.61 


10.17 


16.91 


-.13 


-50.98 


+14.54 


+5.52 


-18.22 


820 


0+8 


8.41 


23 .55 


20.87 


-12.47 


5.88 


33.12 


22.19 


-.13 


-60 .00 


-25.52 


-3.50 


-18.22 


850 


0+8 


8.98 


25.08 


21.76 


-12.78 


6.56 


31.21 


23.35 


-.35 


-58.86 


—16.36 


-4.24 


-16.08 


891 .5 


0+8 


9.71 


23.77 


22.93 


-13.22 


7.25 


26 .40 


24.73 


-.68 


-56.99 


-7.03 


-4.79 


-13.32 


891 .5 





9.71 


11.27 


22.93 


-13.22 


2.20 


9.66 


19.68 


-.68 


-43.55 


+4.29 


+8.65 


-13.32 


900 


P 


9.72 


11.27 


23.04 


-13.32 


2.20 


9.66 


19.78 


-.70 


-43.16 


4.29 


8.69 


-12.91 


950 





9.80 


11.27 


23.65 


-13.84 


2.19 


9 66 


20.30 


-.82 


-41 .04 


4.29 


8.92 


-10.63 


1000 





9.88 


11 27 


24.22 


-14.35 


2.18 


9.66 


20.80 


-.94 


-39.14 


4.29 


9.14 


-8.58 


1050 





9.94 


11.27 


24.78 


-14.83 


2.15 


9.66 


21.27 


-1.04 


-37.41 


4.29 


9.35 


-6.71 


1100 





10.00 


11.27 


25.30 


-15.30 


2.12 


9.66 


21.72 


-1.14 


-35.85 


4.29 


9.55 


-5.00 


1136 





10.04 


11.27 


25.66 


-15.62 


2.09 


9.66 


22.03 


-1.21 


-34.80 


4.29 


9.69 


-3.87 


1150 





10.06 


11.27 


25.80 


-15.74 


2.08 


9.66 


22.15 


-1.23 


-34.42 


4.29 


9.75 


-3.44 


1200 





10.11 


11.27 


26.28 


-16.17 


2.03 


9.66 


22.56 


-1.32 


-33.11 


4.29 


9.93 


-2.01 



Table 13. Thermodynamic functions for Zr Hi. 00 
T in deg K, thermodynamic functions in cal (deg K)-i mole-i. Subscript 1 refers to Zr(a). Subscript 2 refors to H2(g). 































T 


Phases 
present 


H°-H° 

298.15 


c°, 


S° 


G°-H° 

298.15 


Hi-H° 


(Cph 


Si 


Gi-Gf 


Ih-Hi 


(C P h 


S 2 


(? 2 -(? 2 ° 




T 






T 


T 






T 


T 






T 


298.15 


a+8 





7.18 


9.16 


-9.16 


0.02 


6.32 


9.31 


-0.00 


-139.12 


1.71 


-0.29 


-107.62 


300 


a+8 


0.04 


7.20 


9.20 


-9.16 


.02 


6.32 


9.34 


-.00 


-138.29 


1.74 


-.28 


-106.76 


350 


a+8 


1.11 


7.85 


10.36 


-9.25 


.02 


6.42 


10.33 


-.00 


-119.20 


2.86 


+ .07 


-86.95 


400 


a+8 


2.00 


8.70 


11.46 


-9.46 


.02 


6.58 


11.19 


-.00 


-104.73 


4.24 


.54 


-72.02 


450 


a+8 


2.80 


9.60 


12.54 


-9.74 


.03 


6.81 


11.98 


-.00 


-93.32 


5.57 


1.12 


-60 .36 


500 


a+8 


3.52 


10.40 


13.59 


-10.07 


.04 


7.12 


12.71 


-.00 


-84.07 


6.56 


1.76 


-51 .03 


550 


a+8 


4.18 


11.04 


14.62 


-10.44 


.08 


7.51 


13.41 


-.01 


-76.44 


7.05 


2.41 


-43 .38 


600 


a+8 


4.77 


11.49 


15.60 


-10.83 


.13 


7.97 


14.08 


-.02 


-70.06 


7.03 


3.03 


-37.01 


650 


a+8 


5.30 


11.85 


16.53 


-11.23 


.21 


8.48 


14.74 


-.04 


-64.69 


6.75 


3.58 


-31 .62 


700 


a+8 


5.78 


12.39 


17.43 


-11.64 


.31 


9.01 


15.39 


-.06 


-60.09 


6.76 


4.08 


-27 .00 


750 


a+8 


6.26 


13.53 


18.32 


-12.06 


.43 


9.53 


16.03 


-.08 


-56.07 


8.00 


4.58 


-23.00 


800 


a+8 


6.78 


15.93 


19 .26 


-12.48 


.55 


10.00 


16.66 


-.11 


-52.41 


11.86 


5.20 


-19.19 


820 


a+8 


7.02 


17.44 


19.67 


-12.65 


.61 


10.17 


16.91 


-.13 


-50.98 


+ 14.54 


+5.52 


-18.22 


820 


0+8 


7.79 


20.36 


20.44 


-12.65 


5.88 


33.12 


22.19 


-.13 


-60.00 


-25.52 


-3.50 


-18.22 


850 


0+8 


8.29 


23.04 


21.23 


-12.94 


6.56 


31.21 


23 .35 


-.35 


-58.86 


-16.36 


-4.24 


-16.08 


900 


P+5 


9.11 


22.60 


22.55 


-13.44 


7.35 


25.37 


24.97 


-.75 


-56.58 


-5.52 


-4.85 


-12.78 


950 


0+8 


9.77 


20.60 


23.72 


-13.95 


7.75 


20.07 


26.19 


-1.16 


-54.08 


+1.07 


-4.96 


-9.79 


1000 


0+8 


10.27 


19.04 


24 .73 


-14.46 


7.89 


16.68 


27.13 


-1.56 


-51.58 


4.72 


-4.80 


-7.08 


1050 


0+8 


10.66 


17.99 


25.63 


-14.97 


7.88 


14.76 


27.89 


-1.94 


-49.20 


6.47 


-4.52 


-4.62 


1080 .9 


0+8 


10.86 


17.46 


26.15 


-15.28 


7.83 


13.97 


28.31 


-2.17 


-47.80 


6.98 


-4.32 


-3.22 


1080 .9 


8 


10.86 


11.81 


26.15 


-15.28 


2.13 


9.66 


22.62 


-2.17 


-36.43 


4.29 


+7.06 


-3.22 


1100 





10.88 


11.81 


26.36 


-15.48 


2.12 


9.66 


22.79 


-2.21 


-35 .85 


4.29 


7.14 


-2.59 


1136 





10.91 


11.81 


26.74 


-15.83 


2.09 


9.66 


23.10 


-2.28 


-34.80 


4.29 


7.28 


-1.45 


1150 





10.92 


11.81 


26.88 


-15.96 


2.08 


9.66 


23.22 


-2.30 


-34.42 


4.29 


7.33 


-1.03 


1200 





10.96 


11.81 


27.38 


-16.42 


2.03 


9.66 


23.63 


-2.39 


-33.11 


4.29 


7.52 


+0.41 
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Table 14. Thermodynamic functions for ZrHj.25 
T In deg K, thermodynamic functions in cal (deg K)-i mole- 1 . Subscript 1 refers to Zr(a). Subscript 2 refers to Ih{g). 































T 


Phases 


"■ ''m.\5 


C p 


S° 


G°-TI° 

298.15 


Hi-Hf 


(C„)i 


Si 


0i- Gf 


H 2 -H$ 


(C P h 


S 2 


G 2 -G2 




present 


T 






T 


T 






T 


T 






T 


298. 15 


a+8 





7.39 


9.12 


-9.12 


0.02 


6. 32 


9. 31 


-0.00 


-139.12 


1.71 


-0.29 


-107.62 


300 


a+8 


0.05 


7.42 


9.17 


-9.12 


.02 


0.32 


9. 34 


-.00 


-138.29 


1.74 


-.28 


-106.76 


350 


a+8 


1.15 


8.21 


10.37 


-9.22 


.02 


0.42 


L0.33 


-.00 


-119.20 


2. 86 


+.07 


-86.95 


400 


a+8 


2.10 


9.23 


11.53 


-9.34 


.02 


6.58 


11.19 


-.00 


-104.73 


4.24 


.54 


-72.02 


450 


a+8 


2. 95 


10.29 


12.68 


-9.73 


.03 


6.81 


11.98 


-.00 


-93.32 


5.57 


1.12 


-60.36 


500 


a+8 


3.73 


11.23 


13.81 


-10.08 


.04 


7.12 


12.71 


-.00 


-84. 07 


6.56 


1.76 


-51.03 


550 


a+6 


4.45 


11.92 


14.92 


-10.47 


.08 


7.51 


13.41 


-.01 


-76.44 


7.05 


2. 41 


-43.38 


600 


a+8 


5.09 


12.37 


15.98 


-10.89 


.13 


7.97 


14.08 


-.02 


-70.06 


7.03 


3.30 


-37.01 


650 


a+8 


5.66 


12. 70 


16.98 


-11.32 


.21 


8.48 


14.74 


-.04 


-04. 69 


6.75 


3.58 


-31.62 


700 


a+6 


6. 18 


13.24 


17.94 


-11. 70 


.31 


9.01 


15.39 


-.06 


-60. 09 


6. 76 


4.08 


-27.00 


750 


a+8 


6. 69 


14.53 


18.89 


-12.20 


.43 


9.53 


16.03 


-.08 


-56.07 


8.00 


4.58 


-23.00 


791.5 


a+8 


7.15 


16.77 


19.73 


-12.57 


.53 


9.93 


16.55 


-.11 


-53.02 


10.94 


5.08 


-20.05 


791.5 


8 


7.15 


13.12 


19.73 


-12.57 


2.91 


13.12 


18.94 


-.11 


-56. 82 


0.00 


1.26 


-20. 05 


800 


8 


7.22 


13.12 


19.87 


-12.65 


2.94 


13. 12 


L9.08 


-.15 


-56.29 


. 00 


1.26 


-19.44 


820 


8 


7.36 


13.12 


20.19 


-12.83 


3.01 


13.12 


19.40 


-.22 


-55.09 


.00 


1.26 


-18.07 


845.6 


8 


7.53 


13.12 


20.59 


-13.06 


3.09 


13.12 


19.80 


-.32 


-53.64 


.00 


+1.26 


-16.39 


845.6 


P+S 


7.53 


20. 64 


20. 59 


-13.00 


6.49 


31.62 


23.18 


-.32 


-59.05 


-17.56 


-4.16 


-10.39 


850 


0+8 


7.60 


20. 99 


20. 70 


-13.10 


6.56 


31.21 


23. 35 


-.35 


-58.86 


-16.36 


-4.24 


-16.08 


900 


0+8 


8.39 


21.91 


21.94 


-13.55 


7.35 


25. 37 


24. 97 


-.75 


-56.58 


-5.52 


-4. 85 


-12.78 


950 


0+8 


9. 07 


20. 74 


23.10 


-14.02 


7.75 


20. 07 


26.19 


-1.16 


-54.08 


+1.07 


-4.96 


-9. 79 


1000 


0+8 


9.63 


19.63 


24.13 


-14.50 


7.89 


16.68 


27.13 


-1.56 


-51.58 


4. 72 


-4.80 


-7.08 


1050 


0+8 


10.08 


18.80 


25. 07 


-14.99 


7.88 


14. 76 


27.89 


-1.94 


-49. 20 


6.47 


-4.52 


-4.62 


1100 


0+8 


10. 46 


18.03 


25. 93 


-15.40 


7.78 


13.56 


28.55 


-2.31 


-46. 98 


7.15 


-4. 20 


-2.38 


1136 


0+8 


10. 69 


17.45 


26. 50 


-15.81 


7.68 


12.87 


28. 98 


-2. 56 


-45.49 


7.32 


-3.97 


-0.89 


1150 


0+8 


10.77 


17.22 


26.71 


-15.94 


7.64 


12.62 


29.13 


-2.66 


-44.94 


7.35 


-3.88 


-.34 


1200 


0+8 


11.02 


16.39 


27.42 


-16.40 


7.47 


11.76 


29.65 


-2.98 


-43.06 


7.41 


-3.56 


+1.53 



Table 15. 



Ideal-gas thermodynamic functions for normal hydrogen, H 2 {25 percent para, 75 percent ortho) 
T in deg K, thermodynamic functions in cal (deg K)- 1 mole -1 



T 


H°-H°2n.n 
T 


c ° 

Op 


S° 


G°— -H"°298.15 

T 


T 


J/°-H°298.15 

T 


<V 


S° 


G°-n°mAS 
T 


298.15 





6.891 


31 .209 


-31 .209 


845 .6 


4.537 


7.108 


38 .502 


-33.965 


300 


0.043 


6.894 


31 .252 


-31 .209 


850 


4 551 


7.110 


38.539 


-33. OSS 


350 


1 .026 


6 .951 


32.321 


-31 .295 


865 .3 


4 .596 


7.120 


38 .666 


-34 .070 


400 


1.769 


6.975 


33.251 


-31 .482 


891.5 


4.671 


7.130 


38 .879 


-34.208 


450 


2.347 


6.987 


34.072 


-31 .725 


900 


4 .695 


7.142 


38 .946 


-34.251 


500 


2.813 


6.993 


34 .809 


-31 .996 


950 


4.825 


7.180 


39 .334 


-34.509 


550 


3.192 


7.001 


35.477 


-32.285 


956 .6 


4.841 


7.185 


39 .384 


-34.543 


600 


3.511 


7 .009 


36 .085 


-32.574 


1000 


4.944 


7.220 


39 .702 


-34 .768 


650 


3.780 


7.021 


36.647 


-32.807 


1050 


5 .052 


7.261 


40 .055 


-35.003 


700 


4.011 


7 .037 


37.168 


-33.157 


1080 .9 


, r ).ll(i 


7.291 


40 .267 


-35.151 


750 


4.214 


7 .057 


37.653 


-33 .439 


1100 


5.154 


7.309 


40 .395 


-35.241 


791.5 


1 .363 


7.077 


38.034 


-33.671 


1136 


5.223 


7 .343 


40 .631 


-35.408 


800 


4 .392 


7 .081 


38.110 


-33.718 


1150 


5 .249 


7.357 


40 .721 


-35.472 


820 


4.458 


7.092 


38.285 


-33.827 


1200 


5.337 


7.406 


41 .035 


-35.698 
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